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Abstract. We present A3. 5 cm and A6.2 cm radio continuum maps in total and polarized intensity of the barred 
galaxies NGC 1097 (at 2"— 15" resolution) and NGC 1365 (at 9"-25" resolution). A previously unknown radio 
galaxy southwest of NGC 1097 is reported. Apart from a smooth faint envelope and a bright central region, both 
galaxies exhibit radio ridges roughly overlapping with the massive dust lanes in the bar region. The contrast in 
total intensity across the radio ridges is compatible with compression and shear of an isotropic random magnetic 
field, where the gas density compression ratio is approximately equal to 4 and the cosmic ray density is constant 
across the ridges. The contrast in polarized intensity is significantly smaller than that expected from compression 
and shearing of the regular magnetic field; this could be the result of decoupling of the regular field from the 
dense molecular clouds. The regular field in the ridge is probably strong enough to reduce significantly shear in the 
diffuse gas (to which it is coupled) and hence to reduce magnetic field amplification by shearing. This contributes 
to the misalignment of the observed field orientation with respect to the velocity vectors of the dense gas. Our 
observations, for the first time, indicate that magnetic forces can control the flow of the diffuse interstellar gas 
at kiloparsec scales. The total radio intensity reaches its maximum in the circumnuclear starburst regions, where 
the equipartition field strength is about 60 fiG, amongst the strongest fields detected in spiral galaxies so far. 
The regular field in the inner region has a spiral shape with large pitch angle, indicating the action of a dynamo. 
Magnetic stress leads to mass inflow towards the centre, sufficient to feed the active nucleus in NGC 1097. - We 
detected diffuse X-ray emission, possibly forming a halo of hot gas around NGC 1097. 

Key words, galaxies: magnetic fields - galaxies: individual: NGC 1097, NGC 1365 - galaxies: nuclei - galaxies: 
spiral - galaxies: bars - galaxies: structure - ISM: magnetic fields 



1. Introduction 

Radio polarization observations have revealed basic prop- 
erties of interstellar magnetic fields in galaxies of vari- 
ous morphological types (Beck et al. 119961 Beck 2005). 
Until recently magnetic fields in barred galaxies remained 
relatively unexplored, although they can be expected to 
have interesting properties. Strong non-axisymmetric gas 
flows and large-scale shocks will have a major effect on 
interstellar magnetic fields; velocity gradients may en- 
hance the regular magnetic field (Chiba & Lesch 119941 
Otmianowska-Mazur et al. I1997f) and dynamo action 
should also be strongly affected by the presence of a bar 
(Moss et al. I1998|) . According to gas-dynamical simula- 
tions, a shear shock occurs in the bar where the gas 
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streamlines are deflected inwards (Athanassoula 1992a 
1992b; rmer et al. H555j l. To date, none of the numerical 
models includes the magnetic field in a way that allows its 
influence on the shock to be explored, although it is well 
known that magnetized shocks can develop properties dif- 
ferent from those of hydrodynamic shocks. 

Radio observations of barred galaxies can then provide 
useful insight into their interstellar gas dynamics. In order 
to clarify the nonthermal properties of barred galaxies, we 
have performed a radio survey of a sample of barred galax- 
ies (Beck et al. 2002). Among the twenty galaxies studied, 
the prototypical examples NGC 1097 and NGC 1365 are 
most spectacular. Here we discuss in detail and interpret, 
in terms of gas flow and magnetic field models, the non- 
thermal properties of these two galaxies. 



2 



R. Beck et al.: Magnetic fields in barred galaxies. IV. 



The total [polarized + unpolarized] synchrotron ra- 
dio intensity depends on the strength of the total [regular 
+ random] magnetic field in the sky plane and the en- 
ergy density of cosmic ray electrons. Polarized synchrotron 
emission is produced by cosmic ray electrons in the pres- 
ence of either a regular (coherent) magnetic field 1 , or an 
anisotropic random (incoherent) field, or a combination of 
both. Faraday rotation is produced by the component of 
the regular field along the line of sight and thermal elec- 
trons. Since the polarization angle is only sensitive to the 
orientation of the magnetic field in the synchrotron source 
(rather than to its direction), only Faraday rotation is 
sensitive to the distinction between regular and random 
anisotropic magnetic fields. The anisotropy of the ran- 
dom interstellar magnetic field can be significant in barred 
galaxies where it is produced by strong shear and shocks. 
Therefore, the magnetic field strength obtained from po- 
larized intensity, 2 denoted here Bp, must be carefully 
distinguished from that of of the regular magnetic field B. 
Denoting the strength of the anisotropic part of the ran- 
dom magnetic field (b) by b a , we can write Bp = B + 6 2 . 

The first high-resolution radio map of a barred galaxy, 
NGC 1097 (Ondrechen & van der Hulst [Tgg| , showed nar- 
row radio ridges coinciding with the dust lanes, the tracers 
of compression regions along the leading edge of the bar. A 
similar result was obtained for M83 (Dndrechen 11985]) . a 
galaxy with a smaller bar than NGC 1097. The first detec- 
tion of polarized radio emission from a bar was reported 
by Ondrechen (1985) for M83, with a mean fractional po- 
larization of 25%. 

The first high-resolution polarization observations of 
a galaxy with a massive bar, NGC 1097, by Beck et al. 
( 199i5J) had a resolution of 15". Magnetic field orientations 
in and around the bar were shown to approximately follow 
the velocity field of the gas in the corotating frame, while 
the outer field has a spiral pattern. A narrow ridge of 
vanishing polarized intensity indicated deflection of the 
field lines in a shear shock, similar to the deflection of the 
velocity vectors (Athanassoula 1992b). 

Moss et al. 1)2001(1 presented a generic dynamo model, 
based on the model velocity field of Athanassoula (1992b), 
which could explain the major magnetic features observed 
in radio polarization. A similar dynamo model, but now 
based on the specific velocity and gas density fields of 
NGC 1365, is discussed by Moss et al. pHTB) . Beck et al. 
( 2002J presented an atlas of radio maps (in total and po- 
larized intensity) of northern and southern barred galaxies 
observed with the VLA and the ATCA at a resolution of 



1 The regular field is also known as the mean, average, or- 
dered, or large-scale field, but these terms could be misleading 
and are avoided in this paper. 

2 In observational papers it is generally stated that the polar- 
ized emission is a signature of the regular field, without defining 
the term 'regular'. In many papers the possible contribution of 
the anisotropic random field to the polarized emission is ne- 
glected. Our definition is such that only coherent fields are 
called regular. 



30". Harnett et al. I|2004|) discussed ATCA radio observa- 
tions of the peculiar barred galaxy NGC 2442. The circum- 
nuclear rings of the southern barred galaxies NGC 1672 
and NGC 7552 were discussed by Beck et al. (|2005|l . 

In this paper, we present VLA observations of 
NGC 1097 and NGC 1365, with higher resolution than 
in Beck et al. ifT§§§|) and Beck et al. (|2002[) . which re- 
solve the structure of the magnetic field in the ridges of 
NGC 1097 and NGC 1365 and in the circumnuclear star- 
burst regions. Furthermore, we present a ROSAT X-ray 
map of NGC 1097. 

NGC 1097 is a barred galaxy of morphological type 
SBb(s) at about 17 Mpc distance (1" ss 83 pc), with a 
bar of about 20 kpc in length continuing into two optical 
spiral arms. The galaxy plane is inclined by about 45° to 
the line of sight, and its line of nodes has a position angle 
of about —45° with respect to the north-south direction 
(Ondrechen et al. H989|l . with the south-western side closer 
to us. The bar has a similar position angle of about —32° 
and thus lies almost in the plane of the sky. In the galaxy's 
plane, the bar is located at ~ 18° azimuthal angle from the 
major axis. Very little Hi gas has been found in the bar 
(Ondrechen et al. 119891) . CO observations by Gerin et al. 
( 1988J) were not sensitive enough to detect cold gas in the 
bar. Crosthwaite (200Q detected extended CO emission in 
the (1-0) and (2-1) transitions, but the spatial resolution 
was insufficient to see any gas compression in the bar. 
Roussel et al. 12177711 imaged NGC 1097 and NGC 1365 in 
the mid-infrared. Dust emission is relatively weak in the 
bar's ridge, and substantial diffuse emission was detected 
around the bar. 

The nuclear ring of 18" (1.5 kpc) diameter, formed by 
accreting gas, shows enhanced star formation. This con- 
spicuous feature is visible in the optical, CO, radio con- 
tinuum and X-ray spectral ranges. Some spiral dust fila- 
ments were discovered interior to the ring (Barth et al. 
1995] Prieto et al. I2()()5|l . possibly indicating mass inflow 
towards the active nucleus. 

NGC 1097A is a companion dwarf galaxy at about 
18 kpc projected distance towards the north-west (see 
Fig. which causes gravitational distortions in the 
northern half of NGC 1097 and a tidal arm in the Hi 
gas (Ondrechen et al. [1989). 

NGC 1365, at 18.6 Mpc distance (1" « 93 pc), is 
probably the best studied barred galaxy (see the review 
by Lindblad H9"9"9|l . The H I emission and velocity field 
have been studied in detail by Ondrechen & van der Hulst 
lfTMJ|) and Jorsater & van Moorsel fiMfy . The galaxy 
plane is inclined by about 40° to the line of sight and 
its line of nodes has a position angle of about 40° with 
respect to the north-south direction. The north-western 
side is closer to the observer. The bar has a length of 
about 22 kpc and a position angle of about 92°. 

NGC 1365 has an active Seyfert nucleus like 
NGC 1097, but the circumnuclear ring in radio continuum 
is incomplete (Sandqvist et al. I1995f) . Radio continuum 
emission from the bar has first been detected by Beck et al. 
( 2002), it coincides with the dust lanes, as in NGC 1097. 



R. Beck et al.: Magnetic fields in barred galaxies. IV. 



3 



The map of NGC 1365 at A3. 5 cm with 30" resolution 
(Beck et al. 120021 their Fig. 7) shows a deflection of the 
regular magnetic field in its eastern bar, though smoother 
than in NGC 1097. Data at higher resolution and cor- 
rected for Faraday rotation are presented in this paper 
(Sect. EH- 

2. Observations 

2.1. Radio observations 

NGC 1097 and NGC 1365 were observed with the 
Very Large Array (VLA) operated by the NRAO 3 in 
its DnC and CnB arrays at 4.86 GHz (A6.2 cm) and 
8.46 GHz (A3. 5 cm). At both frequencies two IFs, sepa- 
rated by 50 MHz and with a bandwidth of 50 MHz each, 
were recorded. 3C48 and 3C138 were used as primary flux 
calibrators, 3C138 was also used for polarization angle cal- 
ibration; 0240-231 was our phase calibrator. Details of the 
observations are given in Tabled 

Data processing was done with the standard proce- 
dures of AlPS. Visibility data in the same frequency band, 
obtained from the different observation periods, from the 
different arrays and from both IF channels, were combined 
in the uv plane. Maps in Stokes parameters /, Q and U 
were obtained by IMAGR from the combined data at each 
frequency band. Different weightings were applied to ob- 
tain maps with different resolutions (Table [2J ■ Uniform 
weighting (which gives the same weights to each cell in the 
uv plane) reveals the best angular resolution, but larger 
r.m.s. noise than 'natural' weighting (which gives the same 
weights to all antenna pairs and hence emphasizes the in- 
ner uv plane) or 'robust' weighting (an intermediate case). 

At both wavelengths, several values of the 'zero- 
spacing' flux and its weight were tested in the IMAGR 
cleaning process to minimize the negative-bowl effect in 
the total-intensity images caused by missing short-baseline 
visibility data in the uv plane. 

The largest structures visible to the VLA D-array are 
about 3' at A3. 5 cm and 5' at A6.2 cm, less than the full 
extent of the galaxies as observed at A20 cm (Beck et 
al. [2l)()2 ;. No single-dish maps are available to add the 
missing large-scale emission component in total intensity. 
Hence our maps can only be used for spectral index stud- 
ies of the bright ridges and central regions (Figs. and 
We do not expect that any large-scale structure is miss- 
ing in the maps of the Stokes parameters Q and U (and, 
hence, in the polarized intensity maps) because the mag- 
netic field orientation and hence the polarization angle is 
not constant across the galaxies. 

Maps of the linearly polarized intensity (P) and the 
polarization angle (PA) were obtained from the Q and U 
maps. The positive bias in P due to noise was subtracted 
by applying the POLCO correction. 

3 The NRAO is a facility of the National Science 
Foundation operated under cooperative agreement by 
Associated Universities, Inc. 



The final maps of NGC 1097 were smoothed to circular 
FWHM Gaussian beams of 2" (central region, Fig. [11)1 . 4" 
(Figs. H and H2J, 6", 10" and 15" (Figs.HHSJ). Polarization 
angles at 10" resolution were used to compute the map of 
Faraday rotation (Fig. 111( 1 . 

For NGC 1365, no data with the CnB array were ob- 
tained, so that the highest achievable resolution (uniform 
weighting) is 7" at A3. 5 cm and 13" at A6.2 cm. However, 
the maps at these resolutions have low signal-to-noise 
ratios and show little extended emission. In this paper, 
we present only maps obtained with robust and natural 
weighting (see Table 2). Diffuse emission of NGC 1365 is 
best visible in the maps at 15" resolution (Fig. and es- 
pecially in the A6.2 cm maps at 25" resolution (Fig. |3J). 
Maps of the central region at 9" resolution are shown in 
Fig. ^3 Polarization angles at 25" resolution were used to 
compute the map of Faraday rotation fFig. Ill|) . 

The total intensity maps of NGC 1365 were corrected 
for primary beam attenuation. This correction was not 
necessary for the total intensity maps of NGC 1097 and 
the polarization maps of both galaxies because the ex- 
tent of the visible emission is smaller than the diameter 
of the primary beam (5.'4 at A3. 5 cm and 9' at A6.2 cm). 
However, the spectral index is sensitive to small system- 
atic effects, and thus the total intensity of NGC 1097 used 
for the maps of spectral index (Figs. and |§J) had to be 
corrected for primary beam attenuation. 

The values for resolution and r.m.s. noise of the maps 
shown in this paper are given in Table 



2.2. X-ray observations of NGC 1097 

NGC 1097 was observed on Dec 27, 1992, for 2.7 hours 
with the Position Sensitive Proportional Counter (PSPC, 
Pfeffermann et al. 11987(1 on board the ROSAT X-ray satel- 
lite (observation request 600449, PI: H. Arp). We retrieved 
the data from the ROSAT Data Archive and performed a 
standard spatial and spectral analysis with the Extended 
X-ray Scientific Analysis System (EXSAS, Zimmermann 
et al. DM§ . 

To show the extended soft-band X-ray emission, 
the data was spatially binned, exposure corrected and 
smoothed in energy sub-bands (cf. Ehle et al. [1998). 
The smoothing was done with a Gaussian filter with the 
FWHM corresponding to the average resolution of the 
point spread function at the PSPC centre in the individ- 
ual energy sub-bands. The two lowest sub-bands (0.11- 
0.19 keV and 0.20-0.41 keV), smoothed to 52" and 38" 
resolutions, respectively, were added to create a standard 
soft-band (0.1-0.4 keV) map of the X-ray emission of 
NGC 1097 (Fig. EH). 

The spatial distribution and spectral characteristics of 
the soft X-ray emission are discussed in Sect. El 
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Table 1. Observational parameters of VLA observations 



NGC 1097 



NGC 1365 



Frequency (GHz) 
Array 

Pointing (J2000) 
Observing date 
Net observing time (min) 
Shortest baseline 



4.835 & 4.885 
DnC 

02 h 46 m 19!77-30° 16'28'.'55 
1996 May 20 
32 
506A 



4.835 & 4.885 
DnC 

03 h 33 m 35! 50-36° 07'59'.'88 
1996 May 20 and 1999 Feb 19 
20 and 280 
490A and 404A 



Frequency (GHz) 
Array 

Pointing (J2000) 
Observing date 
Net observing time (min) 
Shortest baseline 



4.835 & 4.885 
CnB 

02 h 46 m 19! 87-30° 16' 27'.' 30 
1998 Nov 13 & 15 
560 
765A 



Frequency (GHz) 
Array 

Pointing (J2000) 
Observing date 
Net observing time (min) 
Shortest baseline 



8.435 & 8.485 
DnC 

02 h 46 m 19!81-30° 16'28'.'05 
1996 May 20 and 1997 Oct 07 & 09 
38 and 312 
794A and 704A 



8.435 & 8.485 
DnC 

03 h 33 m 35! 50-36° 07'59'.'88 
1996 May 20 and 1997 Oct 07 & 09 
24 and 311 
790A and 704A 



Frequency (GHz) 
Array 

Pointing (J2000) 
Observing date 
Net observing time (min) 
Shortest baseline 



8.435 & 8.485 
CnB 

02 h 46 m 19!87-30° 16'27"30 
1998 Nov 11 & 14 
560 
1330A 



Table 2. Weights and r.m.s. noise values of total intensity 07 and of polarized intensity (Jq } jj (for Stokes parameters 
Q and U) (in ^Jy/beam) in the final maps 



NGC 1097 NGC 1365 



Beam size 


2" 


4" 


6" 


10" 


15" 


9" 


15" 


25" 










A3. 5 


cm 








Weighting 


uniform 


robust 


natural 


natural 


natural 


robust 


natural 


natural 


07 


9 


7 


8 


11 


17 


12 


12 


20 


a Q,U 


9 


7 


6 


8 


11 


12 


12 


15 










A6.2 


cm 








Weighting 




uniform 


natural 


natural 


natural 




robust 


natural 


07 




12 


8 


15 


20 




14 


14 






12 


8 


9 


13 




14 


14 



3. The radio maps 

3.1. Discovery of a radio galaxy in the field of 
NGC 1097 

A radio galaxy is visible about 90" southwest of the nu- 
cleus of NGC 1097 (Figs. Hand EJ). It is not included in 
the NED and NVSS catalogues. The position of the nu- 
cleus (which is still unresolved at 2" resolution) is RA, 
DEC(J2000) = 02 46 12.75, -30 17 05.3, its flux density 
is 753 ± 10/iJy at A3. 5 cm, and its spectrum is flat (with 
a spectral index a ~ +0.1) which is typical of an active 
galactic nucleus. The northern component is located at 
RA, DEC(J2000) = 02 46 12.7, -30 17 01, its flux den- 
sity is 283 ± 36 /iJy at A3. 5 cm, and its spectum is steep 
(a ~ —0.5). The southern jet shows a spectral steepening 
from a ~ +0.1 to a ~ —1.5 with increasing distance from 



the nucleus. No polarized emission was detected from any 
component. 

3.2. The radio ridges in the bar regions 

The total and polarized radio continuum emission of 
NGC 1097 and NGC 1365 is strongest near the galactic 
centres (see Sects. IET1 and l5~2^l and in radio ridges that ap- 
proximately coincide with the prominent dust lanes seen 
in the optical images (Figs.^®- The lengths of the ridges 
are close to the bar lengths. As gas and stars in the bar 
region rotate faster than the bar pattern in the clockwise 
direction (Ondrechen et al. 1989), the ridges are located on 
the 'downstream' side of the bar which is leading with re- 
spect to the sense of rotation, and the regions 'upstream' 
of the shock front are located on the following side of 



R. Beck et al.: Magnetic fields in barred galaxies. IV. 



5 



Radio Galaxy VLA 6.2cm HPBW=4" 
SO MlcroJy/beam 




02 46 13.5 13.0 12.5 1 2.0 

Rl GHT AS CE NS ION (J2000) 



Fig. 6. Radio galaxy west of NGC 1097 (total intensity), 
observed at A6.2 cm with 4" resolution. 



the bar. The ridges enter the central region tangentially 
to the circumference at the two points where the abso- 
lute maxima of polarized intensity occur, as can be seen 
in the lower panels of Figs, fl] and Pj] Significant polarized 
emission P is also detected outside of the ridges in the up- 
stream regions (Figs. Cj] and and from around the whole 
bar region, forming a smooth envelope (Figs. [3] and C51 see 

The general properties of the radio ridges, summa- 
rized in Table CJl have been determined from the maps 
at A6.2 cm, with 6" (« 500 pc) resolution for NGC 1097 
and with 15" (» 1.4 kpc) resolution for NGC 1365. 

The thermal fractions of the radio emission / t h in the 
ridges were estimated from the observed spectral index 
(Fig. EJ), assuming a constant synchrotron spectral index 
of a s = —1.0. A synchrotron spectral index of a s = — 1 
is expected from particle acceleration in strong shocks 
(a?o = —0.5 intrinsically, steepened by Aa = —0.5 due 
to synchrotron and/or inverse Compton losses). We at- 
tribute the observed flattening of the spectrum with in- 
creasing distance from the centre of NGC 1097 (Table CJJ 
to a larger thermal fraction in the outer parts of the ridge. 
This is in line with the Ha intensity (Quillen et al. H995[) 
which is very small in the inner bar and stronger in the 
middle bar of NGC 1097. In NGC 1365, no similar ef- 
fect was detected. However, the separation of thermal and 



NGC 1097 Spectral Index (3cm/6cm) + Total Int. (3cm) HPBW=10" 
-1.5 -1.0 -0.5 




RIGHT ASCENSION (J2000) 



Fig. 7. Spectral index of NGC 1097 between A3.5 cm 
and A6.2 cm at 10" resolution, shown where the signal 
at each wavelength exceeds 15o7 (to avoid systematic ef- 
fects due to the missing large-scale structures). Contours 
represent A3. 5 cm total intensity at 1,2,3,4, 16,64 times 
150 /^Jy/beam area and the same resolution. 



nonthermal emissions is uncertain because of the possible 
error in spectral index due to the uncertainty in a s , and 
also due to missing large-scale emission in the total inten- 
sity maps. Hence the values for / t h in Table Cj] should be 
regarded as crude estimates. For a more detailed study, 
high-resolution observations at further frequencies or sen- 
sitive, extinction-corrected Ha data are required. 

At the highest available resolution (4", corresponding 
to about 330 pc) the southern total-intensity radio ridge of 
NGC 1097 is resolved into several features (Fig. CJJ which 
are not strictly aligned with the ridge axis, neither do 
they coincide with the dust clouds visible in the optical 
image. The strongest peaks are located upstream of the 
dust lanes, and some of them coincide with star forma- 
tion regions visible in the optical images. According to 
the radio spectral index map (Fig. FjJ) the emission from 
the peaks is a mixture of thermal and nonthermal compo- 
nents. The observed spectral index a reveals strong vari- 
ations along the bar, with values ranging between —0.1 
(expected for purely thermal emission) and less than — 1 
(typical of purely nonthermal emission by cosmic-ray elec- 
trons suffering from synchrotron energy loss in a strong 
magnetic field). 
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The location of the star formation complexes upstream 
of the radio ridges seems to be in conflict with the picture 
of the onset of star formation in the shock, as observed in 
other barred galaxies (Sheth et al. 2002 ). Our observations 
indicate that some compression of gas and magnetic fields 
occurs ahead of the shock front, perhaps in the narrow 
dust filaments mentioned in Sect. 13.31 

The upstream and downstream regions in P are sepa- 
rated by valleys of low polarized emission, the depolariza- 
tion valleys, which are discussed in Sect. 13.31 The polar- 
ized intensities are similar on both sides of the depolar- 
ization valleys, whereas the total intensity is much larger 
on the ridges. Thermal emission cannot account for the 
enhanced total emission in the ridges (see Table |3J). This 
indicates that mainly the turbulent component of the mag- 
netic field is compressed in the shock, while the regular 
field remains almost unaffected. This surprising result will 
be discussed in Sects. |8~T1 and f5~51 

In NGC 1097, the average degree of polarization p of 
the synchrotron emission p — P/I s at A6.2 cm decreases 
from about 30% in the middle southern ridge to about 
5% at the end of the bar (Table |3J), while it is almost 
constant (approximately 20%) along the northern ridge 
(here I s is the synchrotron intensity). In NGC 1365, p is 
low (5-15%) along both ridges. The degrees of polarization 
are generally similar at A3 cm and A6 cm, indicating that 
Faraday depolarization is small (see the right-hand panel 
of Fig. Ill[) . except in the inner ridges and in the central 
regions. In the region upstream of the ridge p is 30-40% 
(in both galaxies), larger than in the ridge. 

The observed half- maximum full width w (in both total 
and polarized intensity - the difference is insignificant) of 
the northern and southern ridges of NGC 1097 increases 
from about 600 pc near the circumnuclcar ring to about 
900 pc near the end of the bar. Assuming a Gaussian ridge 
profile, the observed width w was corrected for smearing 
by the telescope beam (with a Gaussian half-power width 
O) to obtain the intrinsic width wq via Wq — w 2 — 2 . The 
intrinsic width varies from about 400 pc to about 700 pc 
(Table 0). The radio ridges of NGC 1365 are broader and 
shorter than in NGC 1097, with an intrinsic width of ~ 
1000 pc and a length of about 1 '5 (8 kpc). 

The total intensity contrast in NGC 1097 between the 
ridges and the surrounding regions, at ±(7"-12") distance 
from the ridge axis, is about 10 in the inner ridge and de- 
creases to about 5 towards the bar's end (Table 0J). These 
estimates apply to both the northern and the southern 
ridges, although the radio emission in the northern ridge 
is fainter. In NGC 1365, the contrast is about 10 in the 
eastern and 5 in the western ridge, while the total radio 
intensity is higher in the west (Table |3J. 

The contrast in polarized intensity between the ridges 
and their surroundings is quite different from that in total 
intensity (Table 0). St rong polarized emission is observed 
upstream of the radio ridge and dust lane in one half of 
each galaxy (especially to the west of the southern ridge in 
NGC 1097 and south of the eastern ridge in NGC 1365), 
separated from the polarization ridge by a narrow 'depo- 



N SCI 097 3.5cm VLA Total Intensity HPBW=4" 




02 46 24 23 22 21 20 1 9 IS 17 

RIGHT ASCENSION (J2000) 



Fig. 8. Total intensity contours of the central and south- 
ern parts of NGC 1097 at A3. 5 cm at 4" resolution. The 
contour intervals are as in Fig. ^ with the basic contour 
level at 20 /xJy/beam area (the 3cr level). 



larization valley' (see Sect. 13. 3|) . Near the outer end of the 
ridges, the polarized intensity upstream is even stronger 
than in the ridge, so that the contrast is smaller than unity. 
In the northern and inner southern ridges of NGC 1097, 
the polarized intensity is about 7 times higher than in the 
upstream regions. The difference from the total intensity 
ridges suggests that turbulent and regular magnetic fields 
respond differently to compression and shear in the dust 
lane region (see Sect.|HJ). 

The ridge of total intensity observed in NGC 1097 
is systematically shifted from that in polarized intensity 
(Fig. I20|) . As determined from several one-dimensional 
cuts of the southern ridge at 10" resolution, the peak in 
total intensity is about 5" (400 pc) upstream - i.e. to the 
south-west - compared to the peak in polarized intensity. 
Note that the position of the peak of an emission structure 
can be determined with a better accuracy than the beam 
width if the signal-to-noise ratio is high (Harnett et al. 
2004). Figure El shows that the polarized intensity follows 
the optical dust lanes very closely. 

In the northern ridge of NGC 1097, a similar shift be- 
tween the ridges of total and polarized emission is visible, 
though less clearly due to the weaker polarized emission. 
In NGC 1365, an offset has not been detected. However, 
the resolution of our data is lower than for NGC 1097, so 
that a shift of < 10" cannot be excluded. 
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Table 3. Properties of the radio ridges in the bars of NGC 1097 and NGC 1365 from the A6.2 cm data. I is the total 
intensity, P the polarized intensity, ej and ep are the contrasts in I and P between the ridges and the surrounding 
regions, a is the spectral index between A6.2 cm and A3. 5 cm, / t h is the thermal fraction of /, I s = 1(1 — / t h) is 
the synchrotron intensity, and p — P/I s is the fractional polarization of the synchrotron emission. Bj and Bp are 
the equipartition strengths of the total and regular + anisotropic random magnetic fields derived from / s and P (see 
Sect. I3.6(l . The distance along the NGC 1097 ridge is measured from the galaxy's centre. For NGC 1365, average values 
along the ridges are given. 
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Fig. 9. Spectral index of the central region and southern 
ridge of NGC 1097 between A3.5 cm and A6.2 cm at 4" 
resolution, shown where the signal exceeds 5(7/ at each 
wavelength. Contours show the A3. 5 cm total intensity 
at 1,2,3,4,16,64 times 40/xJy/beam area and the same 
resolution. 



3.3. The magnetic field structure in the bar regions 

The B-vectors of polarized emission in the bar region 
change their orientation rapidly upstream of the dust lanes 
in both galaxies (Figs. HB}- This leads to depolariza- 
tion within the telescope beam, observed as a 'depolariza- 
tion valley' parallel to the emitting ridge. This was first 



observed in NGC 1097 by Beck et al. ijTMjjl . Our new 
data show that the valley persists when observed with 
higher angular resolution. Its structure is similar at both 
wavelengths so that it cannot be produced by Faraday 
depolarization. The average distance between the depo- 
larization valley and the ridge in polarized intensity is 
10" (» 800 pc) in NGC 1097 and 11" (» 1.0 kpc) in 
NGC 1365. 

The observed width of the depolarization valley, de- 
fined as the distance between points where the polarized 
intensity drops to half the value outside the valley, is 
about 4" (w 300 pc) in NGC 1097 (Fig. QJ and about 
6" (w 500 pc) in NGC 1365. This is the scale at which the 
regular magnetic field is deflected. Note that a value of the 
width smaller than the resolution is reliable because the 
deflection valley appears only in the map of polarized in- 
tensity (which is not a directly observed quantity). In the 
maps of the observed quantities, the Stokes parameters 
Q and U , the valley is represented by a smooth variation 
with a sign reversal at the location of the depolarization 
valley. The width of the valley depends on the gradients 
in Q and U. Beam smearing decreases these gradients, so 
that the intrinsic width is smaller than that given above. 

According to our A6.2 cm data of NGC 1097 with high 
signal-to-noise ratio (right-hand panel of Fig. EJ, the field 
orientation turns smoothly within a distance of about 20" 
(1.6 kpc) upstream of the ridge. Both depolarization and 
field bending are discussed in Sect. 18.31 In NGC 1365 the 
turning of the field lines is smoother and the depolariza- 
tion valley is broader than in NGC 1097. 

Figure 1101 shows the intrinsic magnetic field orien- 
tations in both galaxies, derived from the £>-vectors at 
A3. 5 cm and A6.2 cm, corrected for Faraday rotation (see 
Sect. 13.4)) and deprojected into the galaxies' planes. 

The intrinsic magnetic field in the southern ridge of 
NGC 1097 follows the dust lane orientation (left-hand 
panel of Fig. I10|l . The orientation of the field and that 
of the dust lane generally agree in the northern bar of 
NGC 1097 as well, but both are less pronounced than in 
the south. As the northern half appears distorted in all 
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spectral ranges, we suspect gravitational distortion by the 
companion galaxy NGC 1097A. In NGC 1365 the field in 
the eastern ridge is along the dust lane, but the resolution 
is insufficient to trace details. 

In the southern bar of NGC 1097 (at 35"-50" or 3- 
4 kpc distance from the centre) the intrinsic pitch angle 
of the magnetic field obtained from the polarized emis- 
sion (i.e., the angle between the field orientation and cir- 
cumference) jumps from about 15° (almost azimuthal) 
in the upstream region (west of the southern ridge) to 
about 75° (almost radial) in the southern ridge (see the 
left-hand panel of Fig. HU|) . which yields a deflection an- 
gle of about 60° . The deflection angle decreases to about 
40° at around 60" radius; this decrease is accompanied 
by a reduction in the contrast in magnetic field strength 
(see Table El • In NGC 1365 (see the right-hand panel of 
Fig. HlUl . the deflection angle is about 70° in the inner 
eastern bar. 

The orientations of the field lines do not everywhere 
follow those of the gas streamlines (in the corotating 
frame) of the hydrodynamic models of Athanassoula 
(|1992b|> and the models for NGC 1365 by Lindblad et al. 
(1996). It is instructive to compare the right-hand panel 
of Fig. EH and Fig. 26b of Lindblad et al. fiMty . where 
it appears that the alignment is reasonably tight in the 
northern part of the galaxy, but not in the southern part. 
This is discussed in detail in Sect. 17.21 

The intrinsic field orientation upstream of the south- 
ern bar of NGC 1097 (left-hand panel in Fig. EH follows 
the small, feather-like dust filaments located upstream of 
the main dust lanes as seen in the optical image (Fig. [31 
around RA, DEC(J2000) = 02 46 21, -30 16 25). The 
reason for this alignment and the origin of the filaments 
are unknown. 

The structure of the incoherent and coherent magnetic 
fields is further discussed in Sects. I£"2.1l and l%"3l 

3.4. Faraday rotation and depolarization in the bar 
region 

The polarization angles at A3. 5 cm and A6.2 cm of 
NGC 1097 and NGC 1365 with 10" and 25" resolution, 
respectively, were used to derive maps of the Faraday rota- 
tion measure RM (measured in rad m~ 2 ) shown in Fig. 1111 
The ambiguity in RM due to the ±7r uncertainty in po- 
larization angle is ±1232 rad m^ 2 . The Faraday rotation 
measure RM is proportional to the product of thermal 
electron density n c (in cm~ 3 ) and the line-of-sight com- 
ponent B\\ of the regular magnetic field (in ^G), integrated 
over the path length L (in pc): RM = 0.81 J n e Bi\ dL. 

RM in NGC 1097 (top left panel of Fig. El is generally 
positive on the western and negative on the eastern side 
of the inner southern ridge, with a jump at the position 
of the depolarization valley, where the orientation of the 
field B± in the sky plane also turns rapidly (see Sect. 13. 3|) . 
This discontinuity in the magnetic field is weaker in the 
middle ridge and vanishes in the outer ridge. However, the 



errors in RM increase with increasing distance from the 
centre due to the decreasing polarized intensities. 

Figure El shows the azimuthal variation of RM in 
two rings in the plane of NGC 1097. In the inner ring 
(left-hand panel), west of the southern ridge (at 190° az- 
imuthal angle), RM reverses its sign and jumps by about 
500 rad m~ 2 , which corresponds to a change of the polar- 
ization angle by about 110° at A6.2 cm and by about 35° 
at A3. 5 cm. At larger radii (right-hand panel of Fig. IT5|) . 
the RM jump decreases to about 200 rad m~ 2 . In the outer 
bar (45"-75" radius, not shown), RM does not jump, but 
changes smoothly by about 300 rad m -2 within a range of 
50° in azimuthal angle. 

Negative RM in the southern bar means that the line- 
of-sight component of the magnetic field points away from 
us. Knowing that the position angle of the major axis of 
the projected plane of NGC 1097 is -45° with the south- 
western side being nearer to us, the radial component of 
the magnetic field in the ridge points inwards (i.e. towards 
the galaxy's centre). 

In the northern ridge of NGC 1097, signal-to-noise ra- 
tios are low, so that RM can be determined only with 
low accuracy. RM values of +100 rad m~ 2 around 330° 
azimuth (Fig. I12fl show that the regular magnetic field 
along the ridge also points towards the centre and has a 
similar strength to that in the southern ridge. 

RM in NGC 1365 (bottom left panel of Fig. [TTJ is 
measured only in the central region fSect. l5~2*)l and a few 
patches between the spiral arms. The polarized emission 
in NGC 1365 is generally weaker than that in NGC 1097. 

Faraday depolarization DP (defined as the ratio of the 
degrees of polarization at A6.2 cm and A3. 5 cm) is shown in 
the right-hand column of Fig. 1111 The amount of depolar- 
ization is small (DP ~ 1) in the outer ridge and upstream 
regions of NGC 1097, but increases towards the centre 
(DP ~ 0.5) where the turbulent field strength is pre- 
sumably larger. Depolarization in NGC 1365 is stronger 
(DP = 0.3-0.6 in the ridges and in the spiral arms) than in 
NGC 1097, possibly due to larger thermal plasma density, 
but the data are of limited quality as polarized intensity 
is low at A3. 5 cm. 

3.5. Field structure in the outer regions and spiral 
arms of NGC 1097 

The optical image of NGC 1097 (Fig.^J shows an elliptical 
structure of about 3f5 x 2' in size, connecting the ends of 
the bar. The ellipse is also the outer edge of the disc- 
like diffuse optical emission. This ellipse can reflect one 
family of orbits in a typical bar potential (see, e.g., Fig. 12 
in Athanassoula 1992a). Highly polarized radio emission 
emerges from along most of this feature (Fig. [3J|, with 
-B-vectors mostly tangent to it. This indicates that the 
regular magnetic field is aligned with the gas flow in this 
region. 

The spiral arms are clearly visible in radio continuum, 
at least the parts which are near to the galaxy's centre. 
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The orientation of the B- vectors at A6.2 cm is along the 
inner northern spiral arm, but almost perpendicular to 
the outer arm in the south-west (Fig.|3J|. Strong Faraday 
rotation seems improbable at such large radii, hence the 
field must be distorted, presumably together with the gas 
flow. 

3.6. The magnetic field strength in the bar regions 

We computed the strengths of the total Bj and regu- 
lar + anisotropic random fields Bp in the radio ridges 
of NGC 1097 and NGC 1365 from the total synchrotron 
intensity / s and the degree of polarization p of the syn- 
chrotron emission, using the maps at A6.2 cm at 10" and 
15" resolution, respectively. (Although at A6.2 cm Faraday 
depolarization may be significant in some regions, we pre- 
ferred to use this wavelength because of the higher signal- 
to-noise ratios than at A3. 5 cm.) I s followed from the to- 
tal intensity by subtracting the thermal fraction estimated 
from the observed spectral index a between A3. 5 cm and 
A6.2 cm (Table 13 . 

We further assumed equipartition between the energy 
of the total magnetic field and that of the cosmic rays (pro- 
tons + electrons), with a number density ratio K of cos- 
mic ray protons to electrons of 100 in the relevant (GeV) 
energy range, a path length through the synchrotron- 
emitting medium of L = 500 pc (similar to the intrinsic 
width of the ridges of NGC 1097, see Sect. EPJl . and a 
synchrotron spectral index of a s = —1. We applied the re- 
vised equipartition formula of Beck & Krause 1(2005(1 . (For 
field strengths beyond 10 /j,G and for a s < — 1 the revised 
formula gives smaller field strengths than the classical es- 
timate.) The results for Bi and Bp are given in Table [31 
These estimates scale with (K/L) 1 /( 3 ~ a *) , so that even an 
uncertainty in K or L of 40% would cause an error in B 
of only 10% (for a s = -1). 

The ridges in total and polarized intensity are pro- 
duced by an increase in Bj and Bp due to compression 
and shear. The strong anisotropy in the random magnetic 
field caused by compression is the main contributor to the 
increased intensity (Sects. HTT1 and I5".2.1(l . 

The observed RM can be used to estimate the strength 
of the component Bii of the regular field B along the 
line of sight if the thermal electron density and the path 
length are known. In principle the thermal electron den- 
sity of the Faraday rotating layer can be derived from 
an estimate of the emission measure of the thermal ra- 
dio emission. The resulting estimates for B» depend on 
the uncertain quantities EM, L and the volume filling fac- 
tor rj of thermal electrons to the power —1/2, and are 
thus more sensitive to the assumed parameters than in 
the equipartition estimate. Bearing this in mind, for an 
estimated EM of 100 cm~ 6 pc in the southern ridge of 
NGC 1097 (where the thermal fraction / t h is about 0.05, 
see Table 01, using Eq. (1) in Ehle & Beck l|1993[) . assum- 
ing that 20% of this emission measure is due to diffuse gas 
(r) = 0.2, as in normal spiral galaxies, e.g. Greenawalt et 



al. 11998(1 . and taking a path length L — 500 pc , we obtain 
n e = (EM77/L) 1 / 2 ~ 0.2cm- 3 . For an RM of lOOradm- 2 
observed in the southern ridge, this gives a strength of the 
regular magnetic field component along the line of sight 
of S|| = RM/(0Mn c L) ~ 1 fiG. 

To compute the full regular magnetic field strength B, 
the angle <f>\ between the field and the sky plane has to 
be known. If the regular field is directed along the ridge, 
the position angle of the field in the sky plane with re- 
spect to the galaxy's major axis is <f>2 — 20°. However, 
the left-hand panel of Fig. ED indicates that the regu- 
lar field may have a position angle of about 45° in the 
galaxy plane, or 02 — 35° in the sky plane. The angle 
<fii between the regular field and the sky plane follows 
from tan <pi = tan <f>2 sin % where i = 45° is the galaxy's 
inclination, and the regular field strength follows from 
B = Bn/shi<f>i. For the range of probable 4>\ values, 
we derive B k, 2 — 4 jiG. This is significantly lower than 
the value obtained from polarized intensity, Bp ~ 12 /iG, 
strongly suggesting that a substantial fraction of the po- 
larized emission is due to anisotropic random magnetic 
fields (Sect. ESDI - 

3. 7. Magnetic field strengths in the spiral arms, disc, 
and envelope 

Substantial radio emission has been detected in the north- 
ern spiral arm of NGC 1097, which emerges from the 
northern end of the bar (around RA, DEC(J2000) = 
02 46 25, -30 14 15, see Figs. □ and EJ. The magnetic field 
orientation suggested by the polarization angles is along 
the arm. The mean degree of polarization is low at A3. 5 cm 
and A6.2 cm (p < 5%). The mean equipartition magnetic 
field strengths obtained from the total synchrotron emis- 
sion (after subtraction of 30% thermal emission) and po- 
larized radio emission in the northern arm are Bi ~ 18 [iG 
and Bp ~ 6[iG, respectively. 

Bright radio emission outside of the bar of NGC 1097 
is also detected in the elongated region emerging from the 
end of the southern bar towards the west (around RA, 
DEC(J2000) = 02 46 20, -30 17 50, see Fig. HJ. Here 
the degree of polarization is quite high (p ~ 10%). The 
equipartition strengths of magnetic fields obtained from 
the total synchrotron and polarized radio emissions are 
Bj ~ 22 /iG and Bp ~ 10 /iG, respectively. This feature 
appears to be a part of the ring-like structure discussed in 
Sect. E21 

The distribution of RM in NGC 1097 does not show 
any systematic variation in the region of the spiral arms 
beyond the bar. However, the RM errors are large due to 
the weak polarized emission. 

The spiral arms of NGC 1365 are brighter than those 
in NGC 1097 both in optical light and in radio continuum 
(Fig. UJ, while the degree of radio polarization is simi- 
larly low (p < 5%) in both galaxies. The mean equipar- 
tition strengths of magnetic fields obtained from the to- 
tal synchrotron and polarized radio intensities are about 
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22 /iG and 3^(G, respectively, for the north-western arm 
of NGC 1365, and about 18 fiG and 5/xG for the south- 
eastern arm. 

In the low-resolution images, both galaxies exhibit an 
envelope of diffuse radio emission (Figs.|^|and[^)l where the 
degree of polarization is high, typically 25% at A6.2 cm. 
Assuming that the emission emerges from an extended 
disc with a path length of 1 kpc, negligible thermal emis- 
sion and a synchrotron spectral index of a s = — 1, we 
obtain Bj ~ 13 ^G and B P ~ 7^G for NGC 1097, and 
Bi ~ 9 AiG and B P ~b^G for NGC 1365. If, however, the 
diffuse emission emerges from a halo with a path length 
of 10 kpc, the fields are weaker by a factor of about 1.8. 

4. The global magnetic structures 

The large-scale structure and direction of the regular mag- 
netic field in a galaxy can be recovered from polarization 
position angles at more than one frequency. We have ap- 
plied a method that seeks to find statistically good fits to 
the polarization angles of synchrotron emission, of a super- 
position of azimuthal magnetic field modes exp(zm<^) with 
integer m, where (p is the azimuthal angle in the galaxy's 
plane. We have used AA3.5,6.2cm polarization angles of 
NGC 1097 (10" resolution) and NGC 1365 (15" resolu- 
tion), averaged in sectors with opening angles in azimuth 
of 10° and 20° , respectively. The nuclear regions, inside 
the innermost rings in Fig. 1131 suffer from strong Faraday 
depolarization (see the right-hand panel of Fig. Ill|l and 
hence represent a different problem for modelling that 
should be addressed elsewhere. 

A three-dimensional model of regular magnetic field is 
fitted to the observations of polarization angles at both 
wavelengths simultaneously. The polarization angle af- 
fected by Faraday rotation is given by ip = ipo + RMA 2 , 
where the intrinsic angle ipo depends on both the regular 
magnetic field and the anisotropic random field (both in 
the sky plane), whereas the second term is only sensitive 
to the regular magnetic field. The vector of the regular 
magnetic field is specified in terms of a Fourier expansion 
in azimuth, and ipo and RM are derived consistently with 
each other. Then the coefficients of the Fourier series are 
fitted to the observed angles using nonlinear least squares 
techniques, and the quality of the fits is verified using sta- 
tistical criteria. The method, and its application to data 
from two spiral galaxies, is described in more detail in 
Berkhuijsen et al. l{TW7|l and Fletcher et al. pUUljl . The 
fitted parameters of the regular magnetic fields in the two 
galaxies are given in Appendix^] The fit parameters given 
there can be used to reconstruct the global magnetic struc- 
tures in the galaxies. The resulting regular magnetic field 
structures of the two barred galaxies are shown in Fig. ^| 
We applied this method to NGC 1097 also in Moss et al. 
(2001), but there used data with lower resolution. With 
our new polarization data we can increase the number of 
sectors, and hence the spatial sensitivity of the model, by 
a factor of 2. 



As discussed in Sect. 13.61 a significant part of the po- 
larized emission can be due to a random anisotropic mag- 
netic field (which does not contribute to Faraday rota- 
tion). Therefore, our fits can be a poor representation 
of the observed polarization angles where the anisotropic 
random field is strong and misaligned with the regular 
magnetic field. The regions where the anisotropy of the 
random field is strong are the radio ridges. The anisotropic 
and regular fields are affected by compression in a differ- 
ent way ('Sect. l8"2~3')) . However, we reasonably expect that 
the anisotropy produced by compression and shear will be 
roughly aligned with the regular magnetic field, as indi- 
cated by comparing Figs. 1131 and 1101 

4.1. NGC 1097 

We were unable to achieve a statistically good fit to the 
data for the inner ring (15" < r < 30", 1.25 < r < 
2.5 kpc) of NGC 1097 using combinations of up to three 
modes selected from m = 0,1,2,3,4 for the horizontal 
magnetic field, combined with either a uniform or a 2tt- 
periodic vertical magnetic field. The reason for this is the 
sharp discontinuity and sign change in the Faraday rota- 
tion measures at the southern end of the bar major axis, 
shown in Fig.ll2land discussed in Sect. 13.41 However, a sta- 
tistically good regular magnetic field model can be readily 
obtained by treating separately the two halves of the ring 
on either side of the bar axis. Then each half can be satis- 
factorily described by the combination of modes m = 0, 2. 
To the north-east of the major axis (azimuth 20° to 190° 
measured from the northern end of the major axis), the 
m = 2 mode is twice as strong as the m = mode. In 
the other half of the galaxy, the two modes have similar 
amplitudes. 

The motivation for splitting this ring into two halves 
is that strong shear in the velocity field in the bar results 
in an abrupt change in sign of the radial component of 
the magnetic field. This type of sudden change, giving 
rise to the discontinuity and sign change in RM, cannot 
be well described by a superposition of a small number of 
azimuthal modes. 

The regular magnetic field pattern required to model 
the AA3.5, 6.2 cm polarization angles - and hence the RM 
discontinuity in Fig. ^]- between 1.25 < r < 2.5 kpc in 
NGC 1097 is shown in Fig. IT31 (left-hand panel). At the 
southern end of the bar, slightly upstream of the bar axis, 
the radial and azimuthal components of the coherent reg- 
ular magnetic field change sign. At the same location the 
pitch angle of the field abruptly increases from pa = 27° 
to pa = 41° . (Note that the change in pa is less than the 
deflection angle discussed in Sect. 13.31 due to averaging of 
the observations in sectors; we lose spatial resolution us- 
ing the model but gain the direction of the regular field.) 
This sharp change in the regular field may be a sign of 
strong shear and compression, and this possibility is dis- 
cussed in Sect. 18.31 Three more reversals in the regular 
magnetic field are also apparent in this ring. Two of these 
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are smoother transitions from positive to negative RM and 
do not produce RM discontinuities (Fig. [T2)l . At the north 
end of the bar there is a stronger and sharper reversal 
that resembles that in the south. However, the weak reg- 
ular field in the north-west of the bar (Fig. ^| means the 
global pattern is less apparent here. 

The second ring (30" < r < 45", 2.5 < r < 3.75 kpc) 
was also modelled by splitting the data into two halves, 
on either side of the bar axis. In this case though, no re- 
versals in the regular magnetic field are required (see the 
left-hand panel of Fig. ll3fl . The variation in RM shown in 
Fig. 1121 fright-hand panel) is not as dramatic as in the in- 
ner ring and the observed sign reversals of RM are consis- 
tent with projection effects due to the galaxy's inclination. 
Our fits for the whole ring are consistent with those for 
the two halves reported here. The regular field in this ring 
is nearly azimuthal, especially in the half of the galaxy 
south east of the major axis. In both halves of the galaxy 
the modelled field consists of a strong axisymmetric com- 
ponent and a weaker higher azimuthal mode. Neither the 
radial nor the azimuthal components of the coherent reg- 
ular field change sign in the model. Although the observed 
RM seem to change sign at <f> ~ 190° (right-hand panel 
of Fig. I12JI . the error bars are large, so that the model 
curve is a statistically good fit. The change in magnetic 
field pitch angle at the southern end of the bar major axis 
is in the same sense as for the inner ring. The change in 
the pitch angle from the purely azimuthal upstream field 
(Pa — 0° ) to the deflected field downstream of the bar is 
~ 40° . 

For the third ring (45" < r < 60", 3.75 < r < 5 kpc, 
we did not split the magnetic field into two halves; there 
is little sign of an RM discontinuity in this ring. Again the 
regular magnetic field is nearly azimuthal in the southern 
half of the galaxy and is deflected after passing through 
the dust lane. The fitted regular field in this ring comprises 
a strong m = mode and two higher modes (m = 1,2) 
with amplitudes which are about half as large. 

4.2. NGC 1365 

Polarized emission in NGC 1365 can be traced to larger 
radii than in NGC 1097, so that more rings can be in- 
cluded in our analysis. The polarization angles in the rings 
shown in Fig. 1131 (right-hand panel) were fitted with a 
combination of m = 0, 1 and 2 modes, except for the third 
ring from the centre where only modes m = 0, 2 were re- 
quired. The fit parameters are given in Table IB~2l There 
is a tendency for pa to become larger (a less azimuthal, 
more radial field) approaching the bar from the upstream 
direction, although this is not the case in all rings. The 
amount of deflection in the field orientation is weaker than 
in NGC 1097. 

A reversal is detected in the regular magnetic field di- 
rection in the north-west quadrant (upper right) of the 
second ring from the centre. This is not seen in the RM 
map (bottom left panel of Fig.lTT|l. but only becomes ap- 



parent when the data is averaged into sectors and mod- 
elled. This reversal is similar to that seen in the ridge of 
NGC 1097. Since the weaker signal from NGC 1365 means 
we had to use only sectors with an opening angle of 20° 
compared to 10° in NGC 1097, it is possible that further 
reversals will be revealed with more sensitive polarization 
observations. 

5. The central regions 

5.1. NGC 1097 

Both galaxies studied in this paper host a central star- 
burst region. The circumnuclear ring of NGC 1097 is a 
prototype example of mass inflow and starburst ignition 
in a bar potential and has been studied in various spectral 
ranges (Gerin et al. 119881 Barth et al.^995 Quillen et al. 
IT335I Perez-Olea & Colina ITMfl Kohno et al. 121^31 . In 
high-resolution infrared images, the ring is resolved into a 
complex network of filamentary spiral features (Prieto et 

The radio continuum ring in NGC 1097 was first stud- 
ied in detail by Hummel et al. ©)87). It appears as an 
almost perfect ring with an average radius of 9"0 (740 pc ) 
(9"5 in north-south and 8'.'5 in east-west directions). 
Although it seems that the ring lies in the sky plane, the 
line-of-sight gas velocity is as large as 230kms _1 at po- 
sition angle —45° (Gerin et al. HHBBI, indicative of the 
ring actually lying in the galaxy's plane and being intrinsi- 
cally elongated along the line of sight. The radio spectrum 
of the ring is nonthermal (see below) so it must contain 
strong magnetic fields. Beck et al. Q19990 detected polar- 
ized emission from this ring. Our new data have much 
higher resolution and also allow us to measure Faraday 
effects. 

The radio map of NGC 1097 at our highest resolu- 
tion, shown in the top left panel of Fig. 1141 exhibits 
several prominent knots in total intensity, mostly coin- 
cident with optically prominent star- forming regions. The 
three brightest knots are separated by almost precisely 
120° in azimuthal angle (being located at position an- 
gles 15° , 135° and 255° , measured anticlockwise from 
the north). No similar periodicity is seen in polarized in- 
tensity. Elmegreen (JT994) proposed that the gravitational 
instability can be responsible for the fragmentation of in- 
ner rings into clouds. 

The spectrum is significantly flatter in the knots of 
the circumnuclear ring (a = —0.65 ± 0.02) than between 
them (a — —0.72 ± 0.02) (see Fig. |§J) which can be due 
to a higher thermal fraction, as suggested by Hummel et 
al. ©87). Ass uming the same synchrotron spectral index 
of a s = — 1 as in the ridges (see Sect. I3.2JI . the thermal 
fraction / t h at A3. 5 cm is 45% in the knots and 35% in 
between. If the observed spectral index is smaller than 
a s , the thermal fraction cannot be determined, and then 
we assume that it is negligible. However, the synchrotron 
spectral index a s in the knots could be larger than — 1, 
e.g. due to a contribution of radio emission from individ- 
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ual young supernova remnants which are expected to have 
an average spectral index of a s ~ —0.6. Furthermore, if 
bremsstrahlung loss of the cosmic-ray electrons dominates 
over synchrotron loss, the synchrotron spectral index is the 
same as the intrinsic one (otQ — —0.5 for strong shocks). In 
this case the observed spectral indices are similar to the 
synchrotron spectral index so that the thermal fractions 
are small. An independent estimate of the thermal emis- 
sion from the optical Ha line is impossible as the quality 
of the existing data is insufficient. 

The equipartition strength of the total magnetic field 
in the knots of NGC 1097 is about 60 /iG (assuming a 
synchrotron spectral index of a s = —1.0 and a path length 
of 500 pc ) , applying the revised equipartition formula by 
Beck & Krause (2005J . The average total field strength 
in the central ring is about 55 /iG. The field strength is 
somewhat larger than that given by Beck et al. (1999) 
because here we assumed a steeper synchrotron spectrum 
and a smaller path length. The degree of polarization in 
the knots is low (1-3%), yielding a regular + anisotropic 
random field of Bp ~ 10 /iG (Table 0J. Alternatively, for 
a s = —0.7 and generally negligible thermal emission, the 
total field strength in the knots increases to about 70 /iG, 
while Bp decreases to 4-8 /iG. 

At the positions where the radio ridges and dust lanes 
enter the central ring of NGC 1097, the degree of polar- 
ization is much higher (note maxima in north-east and 
south-west in Fig. 1141 bottom row panels) which yields 
B P ~ 20 /iG (Table 

The regular (or anisotropic random) magnetic field in 
the central region of NGC 1097 is of a spiral shape and 
extends well inside the ring (see the top right panel of 
Fig. I14|l . The equipartition strength of the magnetic field 
obtained from polarized emission, Bp, is about 15 /iG in- 
side of the ring (north of the centre) (Table HJ. 

Faraday rotation measures are generally small in the 
ring of NGC 1097 (top left panel of Fig.EJ- The RMjump 
in the southern ridge (see Fig. ll2|l continues into the ring. 
With such small RM values, it is impossible to separate 
the coherent and anisotropic random field components as 
attempted in the southern bar region (Sect. 003) ■ We sus- 
pect that, in contrast to the ridges, most of the polarized 
emission from the ring originates in a coherent magnetic 
field, as the conditions for dynamo action are perfect: 
high star-formation rate and strong differential rotation 
velocity of order 330 km s -1 (Gcrin et al. 1988, Kohno et 
al. E003). Note that the average pitch angle of the spiral 
field, corrected for Faraday rotation, is about 39° ± 3° , 
which is larger than in the discs of typical spiral galaxies. 
This fact also supports strong dynamo action. However, 
small RM and strong Faraday depolarization (see below) 
in our present data prevent a search for global modes in 
the structure of the coherent regular field. New observa- 
tions with higher sensitivity at shorter wavelengths (where 
Faraday depolarization is weaker) are required. 

Prieto et al. ( 2005 ) interpret the spiral magnetic field 
in the central region of NGC 1097 as indication of a flow 
of warm gas far away from the galaxy's plane, crossing the 



circumnuclcar ring. Detailed measurements of Faraday ro- 
tation are necessary to clarify the geometry of the mag- 
netic field structure. 

Faraday depolarization between A3. 5 cm and A6.2 cm 
is strong in the ring of NGC 1097 (Fig. ^2 an d compare 
Fig. HU bottom row panels), evidently a result of strong 
turbulent fields and high thermal electron density. Only 
the western tangential point has the same polarized in- 
tensity at both wavelengths. The ridge emerging on the 
western (near) side lies in front of the ring where it cannot 
be depolarized by the medium in the ring. The hot ther- 
mal gas of the halo detected in X-rays is too thin to cause 
significant Faraday depolarization at A6.2 cm (Sect. [BJ . 

The nucleus of NGC 1097 has a flat spectrum with a 
spectral index of a ~ —0.05, typical for active nuclei. An 
investigation of its properties is beyond the scope of this 
paper. 

5.2. NGC 1365 

NGC 1365 has a central starburst region which is even 
brighter in radio continuum than the central ring in 
NGC 1097, and hence the total equipartition fields are 
even stronger. The peak value of about 63 /iG in the inner 
part of the massive dust lane located south of the centre 
is one of the largest field strengths found in any normal or 
barred spiral galaxy so far. Stronger fields were observed 
only in the circumnuclear rings of the southern barred 
galaxies NGC 1672 and NGC 7552 (Beck et al. [MM- Th e 
degree of polarization in the inner dust lane (Table UJ) is 
lower than that in the ring of NGC 1097 (Table EJ). In 
the outer part of the same dust lane, south-west of the 
centre, and in the outer part of the eastern dust lane, the 
degrees of polarization increase significantly. The pitch an- 
gle of the regular field in the ring (radial range 7'.'5-22'.'5, 
or 0.7-2.0 kpc) varies strongly between and —70°. 

RM shows a jump in the central region of NGC 1365, 
south of the nucleus (bottom left panel of Fig. 111(1 . A 
similarity to the field reversal in the ridges of NGC 1097 
and NGC 1365 is possible, but needs investigation with 
future data at higher resolution. 

6. The X-ray halo of NGC 1097 

In order to determine the strength of regular fields from 
Faraday rotation measures, independent information is 
needed about electron density (Sect. ESI, e -g- from ther- 
mal radio or Ha emission which is significant from gas at 
a temperature at a few 1000 K. Hot g observed with 
X-rays, may also contribute to Faraday effects. 

The soft X-ray emission from NGC 1097 observed with 
ROSAT (Fig. 116(1 peaks at the galactic nucleus, which 
could be expected as NGC 1097 is known to harbour an 
active Seyfert nucleus. In addition to the nuclear emis- 
sion (which is smoothed out to a galactocentric radius of 
about 1', corresponding to the point spread function of the 
PSPC at low energies), the soft X-ray emission has con- 
tours elongated along the spiral arms and the bar, and was 
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Table 4. Properties of the central parts of NGC 1097 (r < lkpc) and NGC 1365 (r < 2kpc) from the A3. 5 cm data. 
The locations of the regions are described in the text. The thermal fractions were determined with the assumption of 
a synchrotron spectral index of a s = —1.0. 
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< 0.05 
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0.45 
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Inner dust lane (S) 


0.35 


0.5 


63 
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Outer dust lane (SW) 


< 0.05 


4 


38 
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Outer dust lane (E) 


< 0.05 


22 


21 
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even detected outside the bar. This indicates the presence 
of extended hot gas in a disc and/or in a halo. 

The soft X-ray spur in the north-west points towards 
the companion galaxy NGC 1097A (Ondrechen et al. 
1989J and might be a signature of ongoing tidal inter- 
action between the two galaxies. 

Our successful spectral fit to the dominating X-ray 
emission (background-corrected) from the central region 
of NGC 1097 within a galactocentric radius of 1' used 
a combined emission model including a power law com- 
ponent (to describe the nucleus and unresolved point- 
like sources) and a Raymond-Smith plasma with 'cos- 
mic' metal abundances (Raymond & Smith [T~977[l . both 
affected by Galactic foreground absorption (with hydro- 
gen column density of 1.87 x 10 20 cm -2 ). With the fit- 
ted temperature of 9.6 x 10 6 K and photon index s of 

I. 3 (where the X-ray flux density is proportional to E~ s ), 
the luminosity of the inner part of NGC 1097 in the to- 
tal ROSAT energy band (0.1-2.4 keV) was found to be 

II. 8 x 10 40 ergs -1 , with about 84% of it originating from 
the power law component. 

Unfortunately, the X-ray count rate statistics of the 
extended component outside a radius of 1' (where we ex- 
pect that the emission is no longer affected by the Seyfert 
nucleus) were too poor to allow a separate fit of this emis- 
sion component. Higher sensitivity X-ray observations of 
NGC 1097, e.g., with XMM-Newton, are needed to further 
constrain this unresolved emission. Assuming a Raymond- 
Smith plasma with T = 2.3 x 10 6 K (which seems to be 
typical of the X-ray emitting hot gas in spiral galaxies) 
and correcting for the background and the Galactic fore- 
ground absorption, the count rate of 2.86x 10 -2 counts s -1 
in the outer area corresponds to an X-ray luminosity of 
1.65 x 10 40 ergs -1 in the energy range 0.1-2.4 keV. 

The plasma density n c of the extended component can 
be calculated under the assumption of radiative cooling 
and ionization equilibrium (Nulsen et al. 1984), where 
the luminosity in soft X-rays is given by Lx = 1.11 x 
A(T) n 2 V r/, with A the cooling rate, V the emitting vol- 
ume and r\ the filling factor of the emitting gas. For our 
assumed gas temperature, Raymond et al. (1,976) give a 
cooling coefficient A(T) of about 10 -22 erg cm 3 s -1 . To ad- 



Fig. 16. Soft-band (0.1-0.4 keV) X-ray emission from 
NGC 1097, observed with the ROSAT PSPC detector. 
Contour levels are 2 n / 2 x 3 • 10 4 counts s 1 arcmin 2 
(n = 1,2,...) above the background of 12 • 10 -4 counts 
s -1 arcmin -2 . The background optical image was kindly 
provided by Halton Arp. 

dress the possibility that the emitting plasma is located 
in the disc of NGC 1097 we assumed a cylindrical vol- 
ume (inner radius 1', outer radius 4' and a disc path 
length of 1 kpc) to obtain an electron volume density of 
n c ~ 3.4 x 10 -3 77 -1 / 2 cm -3 . If the thermal energy density 
of the hot gas is equal to that of the magnetic field, as is 
the case in other diffuse phases of the ISM, the total field 
strength should be 5.2 /xG x 77 -1 / 4 . The total field strength 
in the interarm space between the bar and the northern 
spiral arm is estimated from the total radio emission as 
about 13 /iG (assuming equipartition with the cosmic rays 
and an average path length of 1 kpc). Thus, either the fill- 
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ing factor of the hot gas is as small as 3%, or the magnetic 
field energy density is larger than the thermal energy den- 
sity, or the magnetic field obtained from radio emission is 
overestimated. 

Note, however, the asymmetry of the extended soft X- 
ray emission, which is stronger and falls off more slowly 
on the north-eastern (far) side of the galaxy than on 
the south-western (near) side (Fig. 1 1 fc>|) . This may indi- 
cate that a large fraction of the extended X-ray emission 
emerges from a quasi-spherical halo of hot gas in front 
of an absorbing (cool) gas disc. Assuming in this case 
a halo radius of 10 kpc (and excluding again the inner 
radius out to 1'), we obtain an electron volume density 
of n c ~ 2.3 x 10 -3 r\~ x l 2 cm -3 . A total magnetic field 
of 5 /zG x ?y -1 / 4 would be in equipartition with thermal 
energy. If the diffuse radio emission around NGC 1097 
emerges from a halo, the total field strength is about 7 fiG 
(Sect.ETJ), which agrees well with the above estimate from 
the soft X-ray emission for a filling factor T) ~ 0.3. 

As the degree of polarization is about 25% (see 
Sect. I3.7fl . about half of the total field in the halo (or 
disc) is coherent (if the anisotropic random field is neg- 
ligible). The expected Faraday rotation for a coherence 
length of 1 kpc is about 5 x r^ 3 / 4 radm -2 , which is much 
too small to explain the observed Faraday rotation. We 
conclude that the hot halo gas around NGC 1097 does 
not contribute to Faraday effects. 

7. Models of gas flow and magnetic field in the 
bar region 

In this and the following sections, we interpret the radio 
maps in terms of magnetic field models and discuss their 
compatibility with gas-dynamical models of barred galax- 
ies. 

The most important assumption of our interpretation 
is the existence of a shock. Note that shear shocks in bars 
behave differently from classical shocks, and the compres- 
sion region may extend deeply into the downstream region 
(Syer & Naravan ll993|l . No 3-D models of shear shocks in- 
cluding magnetic fields are available yet. 

The best observational tool to identify the location of 
the gas shock is the velocity field of the cold, dense gas. 
In NGC 1097 Hi line emission from the bar is very weak 
so that the velocity field is known only in the outer spi- 
ral arms (Ondrechen et al. 19893 ). CO line emission from 
the bar is also weak, and the resolution is insufficient to 
resolve the shock front (Crosthwaite 2001). Velocity fields 
in the bar region were published for a few other galax- 
ies. Steep velocity gradients occur across the dust lanes in 
NGC 1365 (Jorsater & van Moorsel lTMBI Lindblad et al. 
lilMjl and in NGC 1530 (Reynaud & DownesQ]SH|) which 
indicate that the deflection and compression regions in- 
deed coincide. Sensitive observations of the velocity field 
of NGC 1097 and NGC 1365 in the H I and CO lines are 
required to localize the shock fronts and to measure their 
compression ratio and their rate of deflection in the flow 



direction. Hence there is sufficient evidence for shocks in 
both galaxies of our study. 

Perhaps the most notable features of gas flow in barred 
galaxies are dust lanes extended along the leading edge 
of the bar, which are identified with large-scale shocks 
in the interstellar gas (Athanassoula [1992bJ. Strong gas 
compression and velocity shear are expected to occur in 
the dust lanes. Therefore, it is not surprising that they 
host enhanced magnetic fields and appear as bright radio 
ridges in the radio maps discussed in Sect. 13.21 

Our observations have revealed several unexpected fea- 
tures of both regular and random magnetic fields which 
are similar in both galaxies and may be characteristic of 
barred galaxies in general. The regular magnetic field does 
not seem to be perfectly aligned with gas streamlines in 
the reference frame corotating with the bar. In both galax- 
ies, the regular magnetic field upstream of the bar is ap- 
proximately perpendicular to the dust lanes. The field 
then apparently begins to turn and becomes parallel to 
the dust lanes at their position (see Sect. 13.31 and Figs. [21 
and [SJ) . Hydrodynamic models of the gas flow in generic 
barred galaxies (e.g. Athanassoula |1992b|) have stream- 
lines that behave quite differently, with an acute change 
in the direction of the flow at the leading side of the bar. 
Two possible explanations for this effect are discussed in 
Sect. O and O 

In Sect. 18.11 we show that the enhancement in total 
radio emission observed from the radio ridges is compati- 
ble with compression and shearing of a random magnetic 
field in the shock at the dust lanes. It is striking, how- 
ever, that the observed enhancement in polarized inten- 
sity indicates that the regular magnetic field avoids any 
significant amplification by compression and shear in the 
shock. This behaviour is surprising indeed as the effects 
of compression and shearing are clearly seen in total radio 
emission, as well as in the gas distribution and kinemat- 
ics. In Sect. I8.2~3l we attribute this peculiar behaviour of 
the regular magnetic field to the effects of conversion of 
atomic gas to molecular form in the shock (i.e., the dust 
lane) . 

7.1. The effect of inclination and beam smoothing 

The apparent offset of the depolarized strip from the dust 
lane and the early onset of the magnetic field deflection 
could be a geometrical effect due to the inclination of the 
galaxy's disc and the integration of the emission along 
the line of sight through the disc. Here we demonstrate 
the potential consequences of these effects using synthetic 
maps produced from a simple model of the galaxy in which 
the magnetic field sharply changes its direction by 90° in 
the bar. 

We prescribe a magnetic field in the plane of the model 
galaxy in cylindrical coordinates (r, </>, z), with (j) mea- 
sured in the galaxy's plane from the northern and west- 
ern end of the bar major axis to model NGC 1097 and 
NGC 1365, respectively, and the positive z-direction being 
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Fig. 17. The configuration of the magnetic field used to 
explore the effects of geometry and beam smoothing. The 
bar is located at y = 0. 

towards the observer. The mid-plane horizontal magnetic 
field (B r ,B,p) is shown in Fig. 1171 and is defined as: 



B a 



-B (Z)4>/7T, < <f> < 7T , 

Bq(z){\ - 0/tt), tt < (j) < 2tt , 

B ( Z )(cb/TT- 1), 0<4><7T, 

B o (z)(0/7r-2), 7r<0<27r, 



B 2 = 



The dependence of the modelled field on z is described by 



Bq(z) ex exp(—\z\/hi 



(1) 



with the assumed scale height hg = 4 kpc . 

A z-dcpcndcnt distribution of the form is adopted 
for thermal electron number density n e , with the hs re- 
placed with h c — 1 kpc and the radial distribution trun- 
cated at r = 12 kpc. The number density of cosmic rays 
is assumed to be uniform and the synchrotron spectral in- 
dex is assumed to be a s = — 1. After rotating the major 
axis position angle and inclining the galaxy to the line of 
sight by the same angles as in the galaxies NGC 1097 and 
NGC 1365, we calculate the Stokes parameters Q and U by 
integrating emission along the line of sight over the inter- 
vals ±2/is centred at the galaxy's midplane. We convolve 
the result with a Gaussian of FWHM 10" and 15" in the 
plane of the sky for NGC 1097 and NGC 1365, respec- 
tively. Depolarization occurs due to differential Faraday 
rotation along the line of sight and due to beam smear- 
ing. Figures El an d El show the obtained synthetic po- 
larized intensity in grey scale with apparent polarization 




i i i i T*T"i i I I 

-10 12 

x (arcmin) 

Fig. 18. Synthetic map obtained by rotating the field 
shown in Fig. 1171 with the vertical dependency given by 
Eq. Q to the orientation of NGC 1097. The cosmic ray 
electron distribution is assumed to be uniform and the 
Stokes parameters Q and U are integrated along all lines of 
sight, prior to smoothing with a Gaussian beam of FWHM 
10". The grey scale shows polarized intensity in arbitrary 
units (with darker shades corresponding to larger values), 
and the vectors represent the orientation of the polariza- 
tion plane. The continuous solid line shows the bar major 
axis. 



S-vectors superimposed. Not surprisingly, the abrupt turn 
of magnetic field by 90° results in strips of small polarized 
intensity, similar to the depolarization valleys discussed in 
(Sect. EH ■ 

The combination of inclination and smoothing to the 
beam resolution leads to valleys in P parallel to the bar's 
major axis which are offset by 15" from the ridge of maxi- 
mum P , consistent with the observations (Sect. ET^ . The 
width of the valleys of one Gaussian FWHM is also con- 
sistent with the observations. Because of the inclination of 
the galaxies, the two valleys in each model are not equally 
deep, but the difference is unimportant compared with the 
intrinsic asymmetries found in the observations. The B- 
vectors turn by about 90° in front of the bar, much more 
sharply than in the observations, especially in NGC 1365. 
This implies that the turn of magnetic field near the de- 
polarization valley in real galaxies is smoother than that 
in this model, and so it is partially resolved in our obser- 
vations. 

We conclude that the effects of disc inclination and 
beam smoothing can explain the observed offset of the 
"depolarization valleys" from the ridges, but the observed 
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Fig. 19. Same as Fig.^J but for the inclination and major 
axis orientation of NGC 1365 and a Gaussian beam of 
FWHM 15". 



turning of the B- vectors is much smoother than that pre- 
dicted from our simple model. 



7.2. Misalignment between the velocity field and the 
magnetic field 

As demonstrated above, the smooth deflection of the reg- 
ular magnetic field cannot be the effect of an intrinsically 
sharp field reversal, observed with limited resolution in an 
inclined galaxy, and thus must have a physical reason. 

The velocity field of the dense gas in the bar sharply 
changes its direction in the bar which leads to a shock 
(Fig. I20f) . The smooth turning of the observed I?-vectors 
signifies a misalignment between the velocity field and the 
magnetic field orientation. This may imply that the mag- 
netic field diffusivity is significant. The magnetic field may 
also decouple from the flow of molecular g discussed 
in Sect. 15^31 

The dynamo model of Moss et al. (j2001) predicts sig- 
nificant misalignment (by 20-45°) in the bar region (see 
their Fig. 7). A similar conclusion follows from a dynamo 
model for NGC 1365 (Moss et al. 2005). The misalign- 
ment is a strong indication that the regular magnetic field 
is not frozen into the flow but rather subject to significant 
diffusion. Then the magnetic field can persist on a 10 9 yr 
timcscale only if supported by dynamo action (cf. Moss et 
al.mH). 



8. Magnetic field compression and shearing in the 
radio ridge of NGC 1097 

Here we aim to interpret the observed change in the total 
(/) and polarized (P) radio intensities at the radio ridges 
described in Sect. 13.21 In Table [31 we calculated that the 
contrast in / is ej ~ 5-10 and in P, ep ~ 1-7, where the 
ranges reflect the variation with decreasing distance from 
the galactic centre. First, we estimate how shock compres- 
sion and shearing of the total magnetic field (dominated 
by the turbulent part) is likely to increase /. Then we con- 
sider how compression and shear act to change the regular 
field and hence P. We shall also estimate the contribution 
to P of the anisotropic random field produced from an 
isotropic one by compression and shearing. 



150 



100 



50 



o 
o 

CD 
> 

cc 



-100 



-150 







TV- 'I I, \<- 
\ ft \ 'W' 

\ / 1 1 - X 
I i 1 r \ 


i — ■ — • — • — 




PI6 


/ V 

IP 5 X 

III? \ 








/|: ^ 








I 1 1 

#3' * s / 
/ 3 i ' *" S 
/ T i * 7 ** 


X » . 
_ ^ * \ 


h V 




( b J / 




* A 

\ - " 

fe.- 

■ \ 


U6 J 






\ 

\ 




> ... . 


1 . . . 



1.0 



CD 

0.8^ 



o 

0.6^ 
>. 

0.4.| 

>^ 

0.2 1 

CD 



0.0 



-50 50 

relative azimuthal position [degrees] 



Fig. 20. A cut perpendicular to the bar shock front show- 
ing radial velocity (labelled V r ; thin, black solid), gas den- 
sity (p; thin, black dashed), total (16; thick, grey/red solid) 
and polarized (P16; thick, grey/red dotted) radio intensi- 
ties at A6 cm. The radial velocity and density of the gas 
are from Fig. 10 of Athanassoula (199213 and are derived 
from a generic model of a barred galaxy. The azimuthal 
flow is directed from left to right. The radio data are for 
the southern end of the bar in NGC 1097 and are average 
values for a ring at 30" (2.5 kpc) radius and 15" width. 
The two data sets (model and observations) have been 
aligned using the peak in gas density and the peak (ridge) 
in total radio intensity, defining the zero point of the x- 
axis. Note that this zero point is different from that used 
in Fig. E21 



We treat the region of the dust lanes along the bar's 
leading edge as a shock front where, in addition to com- 
pression, velocity shear results from strong flow towards 
the galaxy's centre behind the shock (Tig. I22|) . We note 
that the offset of the dust lanes towards the leading edge 
of the bar is a sign of a strong shear shock in the veloc- 
ity field (Athanassoula ,1992b). Thus, the magnetic field 
is compressed and sheared at the shock front. Each effect 
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not only amplifies the magnetic field, but also makes its 
turbulent part anisotropic. 

Figure 1201 shows a perpendicular slice through the bar 
region using radial velocity and gas density data from the 
model of Athanassoula l|1992b|l and our radio data from 
NGC 1097. 

We have aligned the two datasets by the peaks in gas 
density and total radio intensity. Figure 1201 illustrates the 
general interpretation of the radio ridges that is developed 
below; compression and shear produce stronger total and 
polarized intensities. The depolarized valley (at —30° rel- 
ative azimuthal position in Fig. I20J) and the strong up- 
stream polarized intensity (at —50°) are also clearly visi- 
ble. 

We assume that the gas speed normal to the shock 
immediately behind it is comparable to the sound speed 
(Roberts et al. 119791 Englmaier & Gerhard[l997) and we 
adopt a sound speed of c s = lOkms -1 ; if the velocity nor- 
mal to the shock is subsonic (e.g. Roberts et al. 0,979), 
the arguments presented below are either strengthened 
or unaffected. In order to simplify our estimates, we as- 
sume that the shock front is in the vertical plane pass- 
ing through the galactic rotation axis (thus neglecting 
the offset of the front from the bar's major axis - e.g., 
Athanassoula [l9"92b|l . so that magnetic field components 
normal and tangent to it are B n = and B t = (B r , B z ), 
where (r, 0, z) are cylindrical polar coordinates in the 
galaxy's frame, centred at the galaxy's centre. 



K = — 2 (see Appendix^) and the gas density contrast of 
an adiabatic shock e p = 4, the model using constant cos- 
mic ray density (Fig. |21(a)| ) predicts e/ ~ 10, as observed 
in the inner bar of NGC 1097 and NGC 1365 (Table EJ. 
Note that, although Fig. [2] is plotted using the geometri- 
cal parameters of NGC 1365, the different inclination and 
bar position angle of NGC 1097 produce similar, but not 
identical, results. The lower observed values of e/ at larger 
radii can be attributed to weaker shear, a weaker shock 
or a combination of both. In order for the equipartition 



model (Fig. 21(b) I to match the observed ej requires both 
a much weaker velocity shear K ~ — 1 and a small in- 
crease in gas density e p ~ 2. Since the velocity and gas 
density fields of these galaxies are not known to a suf- 
ficient accuracy to allow a confident rejection of one of 
the models, we can only note that with a constant cosmic 
ray density, a far greater range of parameters describing 
the shock and shear are compatible with the observed e/, 
whereas equipartition requires both weak shear and low 
compression. 

We conclude that the observed enhancement of total 
radio intensity in the radio ridges is consistent with am- 
plification of the random magnetic field by shock com- 
pression and shear. Better observational data on the gas 
density (CO and Hi) and velocity fields in the bars of 
NGC 1097 and NGC 1365 may provide a novel way in 
which to test the assumption of energy density equiparti- 
tion between cosmic rays and magnetic fields. 



8.1. Contrast in total radio intensity 

In order to estimate the compression ratio (the ratio of 
post- to pre-shock values) of the total synchrotron inten- 
sity, e/, we assume that the contribution of the regular 
magnetic field to the total field is negligible at the re- 
quired level of accuracy, as indicated by our observations. 
In Sect. 13. 61 we estimated the strength of the regular field 
in the ridges as about 2 — 4 ^G while the total field is 
about 25 /iG (Table 0). 

In Appendix ^ we derive two equations for ej when 
an initially isotropic, random magnetic field is subject 
to shock compression and shear. Equation (|A.6J| applies 
when the cosmic ray electron density is not changed by 
the shock (at least, on the scale of our resolution of about 
500 pc) whereas Eq. I|A.7|) includes the common assump- 
tion that the energy density of cosmic rays is equal to 
the energy density of the magnetic field - in this case the 
synchrotron emission is enhanced both by an increase in 
field strength and an increase in cosmic ray electron den- 



sity. Figures 21(a) and 21(b) show how e/ depends on the 
compression ratio of the gas density, e p , for different mag- 
nitudes of the velocity shear, parameterized by the factor 
K. 

The most striking difference between the constant cos- 
mic ray density and the equipartition cases is an order of 
magnitude higher contrast in synchrotron emission pre- 
dicted by the equipartition model. For a shear strength of 



8.2. Contrast in polarized intensity 

Polarized intensity (P) can be increased in the radio ridge 
due to (i) shock compression and shearing of the turbulent 
field producing an anisotropic turbulent magnetic field as 
described above and (ii) compression and shearing of the 
regular magnetic field. These effects will tend to align the 
magnetic field with the shock front and so a strong in- 
crease in P can be expected. 



8.2.1. Polarization from an anisotropic turbulent 
magnetic field 



Equation (|A.8|I gives the degree of polarization produced 
by a compressed and sheared random magnetic field. This 
equation assumes that the cosmic ray density is constant. 
We have not derived a version for the case of energy den- 
sity equipartition between cosmic rays and magnetic field, 
but it is sufficient for the present purposes to note that 
equipartition will tend to increase the degree of polariza- 
tion as discussed in Sect. 5.2 of Sokoloff et al. JT998). 



Figure 21(c) shows the degree of polarization, assum- 
ing po = 0.7 for the intrinsic polarization, for the geom- 
etry of NGC 1365 as a function of e p for different shear 
strengths. Again, the case of NGC 1097 is slightly, but 
not significantly different. In the absence of compression 
(e p = 1), increasing the shear produces a greater degree 
of polarization by producing stronger anisotropy. With in- 
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creasing compression the degree of polarization converges 
to p ~ 0.5, regardless of the shear. However, for the 
strongest shear shown, with K = —8, the behaviour is 
counter-intuitive: stronger compression reduces the degree 
of polarization. This happens because the vertical compo- 
nent of the random magnetic field, b z , is increased by com- 
pression together with b r . The projection of the stronger b z 
into the sky plane has components parallel and orthogonal 
to the sheared 6^, thus reducing the anisotropy of b±. We 
have confirmed this explanation by re-deriving Eq. I|A.8(I 
with b z = 0: in this case, no matter how large \K\ is, the 
slope of p is never negative. 

For e p = 4, the degree of polarization arising from 
anisotropic turbulent magnetic field is p ~ 0.5 in 
NGC 1365 and p ~ 0.4 in NGC 1097, with only a weak 
dependence on the shear strength. These are overesti- 
mates since we do not include any beam depolarization. 
Nevertheless, the expected degree of anisotropy of the tur- 
bulent magnetic field is rather high and should contribute 
significantly to the observed polarized intensity. 

The above estimate is based on the anisotropy in b im- 
mediately behind the shock, as specified in Eq. i|A.l|) . In 
fact the isotropy of the turbulent magnetic field can be re- 
stored rather quickly if the turbulence forcing is isotropic. 
With standard estimates of the turbulent eddy size of 
I = 100 pc and the turbulent velocity v = lOkrns -1 , the 
magnetic field can return to isotropy in the eddy turnover 
time t = 1 1 v — 10 7 yr. Given that, behind the shock, the 
gas speed normal to the shock is equal to the speed of 
sound c s ~ lOkms -1 (Roberts et al. 119791 Englmaier & 
Gerhard 1997] see above) and the width of the radio ridge 
is d ~ 500 pc, the residence time of the gas in the ridge 
is r r ~ d/c s ~ 5 x 10 7 yr, which exceeds the isotropiza- 
tion time severalfold. Hence the degree of anisotropy esti- 
mated above only applies to a region in the close vicinity 
of the shock front and may be lower when averaged over 
our 6" ~ 500 pc telescope beam. Thus, whilst shock com- 
pression of the turbulent magnetic field will give rise to 
a strong increase in / (see Sect. I8.1fl it may only have a 
weaker effect on P. 

8.2.2. Polarization from a sheared regular magnetic 
field 

Now we consider the effect of shear in the velocity field 
on the regular (coherent) magnetic field. Gas-dynamical 
simulations (e.g., Fig. 10 of Athanassoula 1992b, see also 
Fig. E3 here) confirm that the radial velocity in the dust 
lane region changes very rapidly near the ridge, from an 
outward to an inward direction, and subsequently varies 
slowly. If the radial velocity changes by AV r over the dust 
lane of a width d, the radial magnetic field B r produced 
from the azimuthal one by the shear is 

_ _ AV r — AV r — , ^ 

Br^Bf— ^r r ~B -~1OB , (2) 

el c s 

where r r ~ d/c s is the residence time in the shock region, 
AV r — 100 km s -1 and c s ~ lOkms -1 . The enhancement 



of the magnetic field tangent to the shock front by com- 
pression is €b — 4 for an adiabatic shock. This is a factor 
of more than two smaller than the amplification by shear, 
Eq. J2J- I n addition, the regular magnetic field is almost 
perpendicular to the shock upstream of it, which reduces 
the effect of compression and enhances the effect of shear. 
So, we neglect the compression of the regular magnetic 
field. 

A similar estimate of regular field amplification in the 
ridge follows from the observed deflection of the polariza- 
tion plane discussed in Sect. 13.31 where we note that mag- 
netic pitch angles upstream and downstream of the shock 
front in NGC 1097 are pai — 15° and PA2 — 75°, respec- 
tively. Given that B r jB^ = tanpA and B^i w B^2, since 
the azimuthal field component is nearly normal to the 
shock front, we obtain B r 2/B r i ~ t&np^/ tanpAi — 14, 
which is in fair agreement with Eq. (J5J. 

Such an enhancement of the regular magnetic field 
must result in a very significant increase in polarized in- 
tensity. The enhancement factor of polarized emissivity 
can be estimated using Eqs. (|A.4|I and (|A.5|I . but now 
written for the components of the transverse (i.e. in the 
plane of the sky) regular magnetic field B± = (B x /,B y i). 
Assuming for the sake of simplicity that, in the galaxy's 
plane, the regular magnetic field is purely azimuthal in 
front of the shock, B\ = (0, 5^,0), and purely radial be- 
hind it, £?2 = (B r ,0,0), and that the enhancement factor 
in polarized synchrotron emissivity £p 2 /epi is equal to 

2 

that in B ± (i.e., the energy density of cosmic rays is con- 
stant across the shock), we obtain 

e P = 

f AV r \ 2 1 - cos 2 4> sin 2 i 
\ c s J 1 - sin 2 <j) sin 2 i 

f 60 in NGC 1097 , 
~ \ 90 in NGC 1365 . 

If the region with strong shear is not resolved then the 
increase in polarized emission can be smaller than this. 
However, in order to achieve ep j0 bs — 3 (Table EJ would 
require the width of the sheared region to be about d ~ 
We P<oha /e P ~ 50 pc in NGC 1365 (where W = 1.5 kpc is 
the beamwidth) and even less in NGC 1097. In addition, 
the observed e p does not increase in our maps with higher 
resolution, so the low observed values are not due to our 
limited resolution. 

We conclude that the observed contrast in polarized 
intensity is inconsistent with the expected amount of mag- 
netic field enhancement if it is subject to a full amount of 
shearing in the dust lanes. 

8.2.3. Decoupling of the regular magnetic field from 
the molecular clouds 

From the above arguments, it is clear that the com- 
bined effects of anisotropic turbulent magnetic field and 
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sheared regular magnetic field should produce a signifi- 
cant increase in the polarized intensity around the shock 
front/dust lane. However, we observe only a modest con- 
trast of e P ~ 0.5-7 (Table EJ. 

The simplest way to resolve the conflict between the 
expected and observed ep, without having to abandon our 
successful description of ej in Section l^Tl is to suggest that 
the regular magnetic field resists shearing by the radial 
velocity field. As we argue in what follows, this can hap- 
pen if the regular magnetic field becomes decoupled from 
the dense gas (molecular clouds) in the shocked region. A 
typical time scale for the formation of H2 molecules is of 
order 10 6 yr ( Jura 119751 Bergin et al. I2004f> . much less 
than the residence time in the ridge r r ~ 5 x 10 7 yr. As 
molecular clouds form, it is plausible that they become 
detached from the regular magnetic field, carrying with 
them and amplifying only the random, small-scale mag- 
netic field. One mechanism (Ohmic diffusion) has already 
been proposed for this process by Mestel & Strittmatter 
(1967) and others are plausible, such as ambipolar diffu- 
sion and reconnection of the external field as a forming 
cloud rotates; further discussion is beyond the scope of 
this paper (see Fletcher & Shukurov l2005f) . 

In order for the regular field - and hence the polarized 
intensity - not to be increased as strongly as suggested 
by Eq. J2J, it must be sufficiently strong to prevent the 
diffuse gas, to which it is coupled, from being sheared. 
Consider a shear flow V = (V x (y), 0, 0) with a horizontal 
magnetic field embedded into it, B = (B x , B y ,0), where 
the x-axis is parallel to the radial direction and the y-axis 
is directed along azimuth. The magnetic braking of the 
radial shear flow is controlled by the x-component of the 
Navier-Stokes equation, 



dV x 



J_ -g dBx 

dt ~ 4np v dy 



(4) 



where we have neglected the x-derivative of magnetic field 
in comparison with its y-derivative. For simplicity, con- 
sider the effect on the magnetic field of compression alone. 
This will result in a very conservative estimate of the effect 
of magnetic field on the flow since additional amplification 
by the shear is neglected. Thus, B x is compressed over a 
length d\ by a factor e p , whereas B y remains unchanged: 



dB x /dy 



X)B x jd\. Since the azimuthal speed be- 



hind the shock is close to the speed of sound c s (see above), 
the increment in the radial velocity produced by magnetic 
stress over a length c?2 follows as 



AV X 



£ p — 1 B x B y d2 



70 



c s 47rp di 
km / B x B y 



(5) 



s \10/j,G 2 J V0.2cm- 3 / VlOkms- 1 , 

where the numerical value has been obtained for d± = di 
and € p = 4. We have the following additional constraints: 

n c B\\ ~ 0.3/iGcm~ 3 , 

n\L ~ 100rycm" 6 pc, (6) 

n c = Xn, 



where the first is derived from Faraday rotation values 
in the ridge (Sect. 13.4(1 . the second from the estimated 
EM with i] the filling factor of the diffuse ionised gas 
(Sect. I3~B|) . and the third defines the degree of ionisation 
X of the diffuse gas. For the inclination of NGC 1097 and 
NGC 1365, L ~ 2\/2h where h is the scale height of the 
diffuse gas. Combining Eq. JSJ, with B x = B y = B/y2, 
AV X — 100 km s -1 and c s = lOkms -1 , and Eqs. © we 
obtain expressions for the regular magnetic field strength, 
diffuse gas density and the scale height of the ionized layer 
in terms of the least well known parameters 77 and X: 



B 
TjiG 



1 cm 



1 kpc 



4X- 1 /3 ; 

0.iA -2/3 , 
SrjX- 2/3 . 



(7) 



In other words, Eq. (JSJ) - and hence the suppression of 
shear in the diffuse gas - is consistent with the observed 
RM and EM and the assumed disc thickness if Eqs. J7J 
give reasonable values of B, n and h, using the chosen 
fractional ionisation X and the diffuse gas filling factor 
rj. This is the case for X > 0.5 and 77 < 0.2, plausible 
values for both parameters. For example, X — 0.5 and 
rj = 0.2 gives B ~ 5 fiG, n ~ 0.2 cm~ 3 and h ~ lkpc. 
The latter is similar to the scale height of the diffuse gas 
in the Milky Way (Reynolds 1990) and the filling factor 
rj = 0.2 is similar to that of the diffuse thermal electrons 
in normal spiral galaxies (Greenawalt et al. P-998). 

To summarise this section: the increase in polarized in- 
tensity in the ridge, with respect to the upstream region, 
is much lower than expected because the regular magnetic 
field decouples from the dense molecular gas clouds and is 
sufficiently strong to prevent shearing and compression in 
the diffuse ionised gas. The regular magnetic field of about 
5 fiG in strength can be dynamically important in the dif- 
fuse gas of density ss 0.2 cm -3 . If the ionization degree of 
the diffuse gas is of order 50% and the scale height of the 
diffuse ionised gas is about 1 kpc , the resulting emission 
measure is compatible with that observed in the southern 
ridge of NGC 1097. 

8.3. Variation in RM across the dust lane 

The variation of the Faraday rotation measure across the 
inner southern radio ridge region of NGC 1097 (left-hand 
panel of Fig. is unusual. The discontinuous change in 
sign of rotation measure arises because the radial com- 
ponent of the regular magnetic field changes direction 
across the shock front. This counter-intuitive arrangement 
of magnetic field directions is confirmed by the fitting 
of the large-scale polarization structure with azimuthal 
Fourier modes presented in Sect. Figure l2~2l illustrates 
how such a magnetic field configuration can be produced 
by a shear shock. We should remember that the observed 
polarization vectors are dominated by the anisotropic field 
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Fig. 22. An illustration of the shearing of the regular mag- 
netic field in the radio ridges. The orientation of the figure 
and the field direction roughly correspond to those in the 
southern part of NGC 1097. A magnetic line shown with 
continuous line is sheared by a velocity field component 
along the radio ridge (dust lane) region, whose vectors are 
indicated with dotted lines; the thin straight line shows 
the orientation of the original magnetic field; the unre- 
solved part of the system is shown shaded. 



which is strictly aligned along the shock front. The sheared 
regular field sketched in Fig. visible only in Faraday 
rotation, is tilted with respect to the shock front. 

The sharp change in magnetic field orientation at the 
sheared region (Fig. 1221) explains the origin of the depo- 
larized valley: If the region of strongest shearing is narrow 
compared to the beam (as we argued in Sect. I8.2T2*|) . then 
the almost orthogonal field orientations will give rise to 
strong, wavelength-independent beam depolarization. A 
shift of the polarization ridge downstream with respect 
to the total intensity ridge will also result, as observed 
(Sect.E2>. 



9. Mass inflow through the nuclear ring of 
NGC 1097 

The nuclear ring of NGC 1097 is a site of intense star for- 
mation with a rate of 5M Q yr" 1 (Hummel et al.^87) and, 
presumably, strong gas inflow to fuel the central activity. 
Beck et al. (1999) suggested that the stress required to 
drive the inflow into the central region can be provided 
by the magnetic field. The total magnetic field in the nu- 
clear ring and its regular + anisotropic random part are 
estimated, from equipartition arguments, to be Bj ~ 50- 
60fiG and B P ~ 10-20 ^G, respectively (Table©. The 
regular magnetic field has the form of a trailing spiral 
with a pitch angle of pa — 40°. This field produces the 
r0-component of the stress tensor T r $ = —UA r UA4>, where 
ma = -B/(47T j o) 1 / 2 is the Alfven velocity based on the total 
magnetic field. 

The accretion rate driven by this stress is M = 
2-K<j£l~ 1 T r< j ) , where a ~ 2hp is the gas surface density, 



with 2h the disc thickness, and £1 is the angular velocity 
of the gas (e.g., Balbus & Hawley[1998). This yields 

M ~ ~(b^+B r B^) , (8) 

where b and B are the turbulent and regular magnetic 
field components, respectively, and overbar denotes aver- 
aging - cf. Eq. (TJJ . 

With c = brb^/b 2 , we obtain 

M ~ -^(c¥+±B 2 sm2p) 
11 A 

~ — c Mq yr -1 , 

where c ~ -1 is a numerical factor (see below) and 
the estimate has been obtained for h = 0.1 kpc, £1 = 

450kms- 1 kpc" 1 ,62 1/2 = 50 fiG, B = 10 fiG and p A = 
40°; we have retained only the (dominant) first term in 
Eq. (JHJ. Hence, we derive a mass inflow rate of a few 
Mq yr -1 which is sufficient to feed the active nucleus. 

A correlation between the azimuthal and radial com- 
ponents of the turbulent magnetic field used above arises 
from the shearing of turbulent magnetic field by differ- 
ential rotation. From arguments similar to those used to 
obtain Eq. (|A.3() and assuming that the relevant turbulent 
timescale is h/c s (as in accretion discs) we obtain 

d£l h q£lh 
&0 = r— b r ~ b r , (9) 

&T C s C s 

where the rotation curve is assumed to be of the form 
Q, oc r~ q , with q = 1 corresponding to a flat rotation 
curve. Assuming that b 2 . = \b 2 1 we obtain c = b r b^jb 2 ~ 
~^qflh/c s . 

The problem of mass transport from a circumnuclcar 
ring into an active nucleus has been addressed by many 
authors (see review by Knapen fj004 ) . The inflow rate 
due to magnetic stress is much larger than that nor- 
mally produced by graviational torques in the nuclear ring 
(Maciejewski et al. I2002fl . We note that turbulent stress 
can result in a similar inflow rate given that magnetic field 
and turbulence are in energy equipartition. Strong mag- 
netic fields are known to exist in the centres of most galax- 
ies from the observation of strong synchrotron emission. 
Hence magnetically driven inflow can be important not 
only in barred galaxies, but also in normal spiral galaxies 
with a magnetized inner disc (Moss et al. I2000|l . 

10. Summary and conclusions 

— We have detected polarized emission in the galaxies 
NGC 1097 and NGC 1365 which gives important in- 
formation about the interaction of the gas flow and the 
magnetic field structure in barred galaxies. 

— The total magnetic field is strong in the central star- 
forming regions of NGC 1097 and NGC 1365 (about 
60 fiG) and in the radio ridges along the galaxies' bars 
(20-30 /xG). 
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— In both galaxies, the magnetic field orientation changes 
strongly, but smoothly in front of the radio ridges. This 
causes a valley of depolarization where the polarized 
emission almost vanishes. The smooth deflection of the 
magnetic field is different from that of the (likely) ve- 
locity field of the dense gas which changes its direction 
abruptly. 

— The enhancement of total radio intensity in the ra- 
dio ridges along the bars of NGC 1097 and NGC 1365 
can be explained by shock compression and shearing of 
the isotropic random field, together with compression 
of the gas, by a factor of about 4 and constant cos- 
mic ray density across the shock. Local equipartition 
conditions, where the cosmic rays are compressed to- 
gether with the magnetic field, are hardly compatible 
with our data as they require both low compression 
and low shear. 

— The enhancement of polarized intensity is much lower 
than that of the total synchrotron intensity. The regu- 
lar magnetic field resists shearing because it is decou- 
pled from the dense gas, and the Maxwell stress is suf- 
ficient to oppose efficiently the development of a shear 
flow developing in the diffuse gas to which the regular 
magnetic field remains coupled. This also contributes 
to the misalignment between the magnetic field and 
the velocity field of the dense gas in front of the ridges. 

— The enhanced polarized emission in the ridges can be 
explained by the presence of anisotropic random fields 
due to compression of isotropic random fields in the 
shock. These anisotropic fields dominate over the reg- 
ular (coherent) fields which are visible through their 
Faraday effects. 

— Faraday rotation measures reveal jumps across the 
ridge, where the amplitude decreases with increasing 
distance from the galaxy's centre. This behaviour re- 
sults from shearing of the regular field. 

— We have demonstrated that radio continuum data of 
the total and polarized emission form a new and power- 
ful tool to determine the shock strength and the shear 
rate in the shearing shocks of barred galaxies. If the 
gas compression is known from independent measure- 
ments (Hi and CO), the validity of energy equiparti- 
tion between magnetic fields and cosmic rays can also 
be tested. 

— Enhanced polarized intensity ahead of the shock fronts 
of NGC 1097 and NGC 1365, at several kpc distance 
from the bar, could be a signature of regular fields 
generated by dynamo action. Alternatively, compres- 
sion of gas and magnetic fields may occur in the region 
of narrow dust filaments visible on optical images. The 
magnetic field is aligned with these filaments. 

— Magnetic braking in the circumnuclear regions of 
galaxies is a fundamental process which may solve the 
long-standing problem of how to feed active galactic 
nuclei. We have shown that the magnetic fields in the 
circumnuclear ring of NGC 1097 are sufficiently strong 
(~ 60 fiG) to drive a mass inflow of several M yr _1 . 
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- X-ray emission in NGC 1097, detected with ROSAT, 
reveals hot gas (of a temperature about 10 7 K) in the 
central region and also an extended hot gas compo- 
nent. A crude estimate of the electron density in the 
extended component yields n c ~ 3.4 x lO -3 ?? -1 / 2 cm~ 3 
(where r\ is the volume filling factor of the hot gas). 
If the extended X-ray and the extended radio com- 
ponents both emerge from a halo and r\ ~ 0.3, the 
derived thermal energy density of the hot gas is com- 
parable to the energy density of the magnetic field. 
A more detailed spectral and spatial X-ray analysis to 
study further the hot gas in the ISM of barred galaxies 
such as NGC 1097 requires deep observations at high 
throughput observatories like XMM-Newton. 
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Appendix A: The effect of compression and shear 
on a random magnetic field 

Here we derive equations describing the increase in to- 
tal synchrotron emission and polarization resulting from 
shock compression and shearing of an initially isotropic, 
random magnetic field. For the total emission we consider 
two cases: first, and simpler, we assume that the energy 
density of cosmic ray electrons is constant (i.e. the den- 
sity of cosmic rays is not affected by the shock); then we 
present the more complex model where the cosmic rays are 
in energy equipartition with the magnetic field. Finally, 
we examine how much polarized emission emerges from 
the anisotropic field produced by shock compression and 
shearing. 

The component of the random magnetic field tangent 
to the shock front is compressed together with the gas, 
whereas the normal component remains unchanged, 

bt2 = Cpfetl , &n2 = Kl , (A.l) 

where e p — pij P\ is the density ratio and subscripts 1 
and 2 refer to pre- and post-shock values, respectively. 
We assume for simplicity that the random field is isotropic 
upstream of the shock, 

<&> = *<6?i> = !<&?> , (A.2) 

where angular brackets denote averaging. Thus, the com- 
pressed field is anisotropic, with the dominant components 
tangent to the shock front. 

The anisotropy of the random magnetic field in the 
postshock region is further enhanced by velocity shear. 
We introduce cylindrical polar coordinates (r, </>, z) and 
assume for simplicity that the shock is parallel to the r 
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direction. If the radial velocity changes by AV r across the 
shock region of a width d, the radial magnetic field b r 2 
produced from the azimuthal b^i by the shear is 



AV r 
0<hl —T T 



AV r I 



001" 



-2b, 



(A.3) 



where r r = min(d/c s , 1/vq) is the relevant amplification 
time (the minimum of the residence time within the 
sheared region, d/c s (where c s is the sound speed), and 
the eddy turnover time l/v). Since the turbulent veloc- 
ity vq ~ c s and d > I, we obtain the second equality in 
Eq. {0]| , and the th ird follows for AV r = -100 km s^ 1 
(Athanassoula 1992b), vq = lOkms -1 , d = 500pc and 
I = 100 pc. 

The radial component of the random magnetic field be- 
hind the shock thus consists of two parts, the compressed 
one, denoted by br as given in Eqs (|A.1(1 and ljA.21) . 
and the sheared part b r given by Eq. lfO|) . Although 
(brbcf,) ahead of the shock (if the field anisotropy 
can be neglected there), the shearing produces correlated 
field components behind the shock, {b r b^) — {bi^b^) ~ 
(b 2 ,)(AV r l)/(v d) (where AV r < 0). For the post-shock 

field produced by compression, we have = £pb r i and 
likewise for the z-component, whereas the ^-component 
remains unchanged as it is normal to the shock front. The 
compression enhances an isotropic random magnetic field 

1 /2 

by a factor of [|(1 + 2e 2 )] . For a strong shock (e p = 4) 
the compression factor is approximately 3. 

Since synchrotron intensity depends on the magnetic 
field component in the plane of the sky, its compression 
factor £] depends on the orientation of the shock front 
with respect to the observer. Appropriate expressions for 
the magnetic field projection to the plane of the sky, 
b± = (b x i,b y i), in terms of its cylindrical components 
b = (b r ,b,f,,b z ) in the galaxy's reference frame can be 
found in Appendix A of Berkhuijsen et al. (1997): 



b x < = b r cos 4> — b^ sin 4> , 

by/ — (b r sin <fi + b^ cos <f>) cos i + b z sin i , 



(A.4) 
(A.5) 



where i is the inclination angle of the galactic disc (i = 
corresponds to the face-on view) and <f> is the azimuthal 
angle in the galaxy's plane measured anticlockwise from 
the northern end of the major axis, related to that in the 
sky plane, 0', by tan <f> cos i = tan<^'. The x'-axis in the 
sky plane points to the northern end of the galaxy's major 
axis; here and below, coordinates in the plane of the sky 
are denoted with primes, (x',y'). 

In order to calculate the contrast in synchrotron inten- 
sity, we consider the two models described below. 

A.l. Total intensity contrast for uniform cosmic ray 
distribution 

If the energy density of cosmic rays does not vary across 
the radio ridge (i.e., at a scale of 1 kpc), the increase 
in synchrotron emission depends solely on the change in 



6j_, the magnetic field component parallel to the plane of 
the sky, and any variation in the disc thickness between 
the shocked and upstream region. Then Eqs. (|A~T|l . ljA~2)l 
and l|A.3j) yield the contrast in total synchrotron intensity 
between the shocked and upstream regions: 



(b 



2 ) 

±21 



(&ii> In 
= I ^[2+^ + ^-1X1 

-K(K + sin 2<f>) sin 2 i] , 



■ 2 ■ 2 

sin i cos 



<t>) (A.6) 



where K = (AV r l)/(vod) is a measure of the strength 
of the shear and for the canonical values given above 
K ~ — 2. We have assumed that different components of 
the compressed part of b remain statistically independent, 
(br^b^) — and similarly for the other components. For 
Li = L 2 , e p = 1 and AV r — 0, we obtain, as it should 
be, ex = 1; for i = (the face-on view), the contrast in / 
becomes independent of azimuth. 



A. 2. Total intensity contrast under energy equipartition 
between cosmic rays and magnetic fields 

Here we consider the case where the energy density of cos- 
mic rays is proportional to the energy density of the total 
magnetic field. Then the synchrotron emissivity behind 
the shock will be increased both because of the field am- 
plification by compression and shear and, additonally, due 
to the accompanying increase in the cosmic ray electron 
density where n CI oc b 2 . Now we obtain, for L\ = L^: 



K 2 b 2 ) 
(b 2 ±1 b 2 ) 

1 + ^(e 2 + K 2 - l) 2 (cos 2 i + cos 2 <?(> sin 2 i) 

„2 „■ , n 2 i -2 



^(e 2 + K 1 - 1)(4 + 9 cos^ i + 7cos^ 4> sin^ i) 



+ ~ !)( e P + R2 - !)(! + cos2 * sin2 
+ ±(e 2 p -l)(3e 2 p + 4) sin 2 i 



e z p + K z ) sin 2(j> sin 2 i , 



(A.7) 



following careful treatment of the averages of the compo- 
nents of the field and their cross terms ((&*), (tftb^) etc.) 
and some algebra. 

A.3. Polarization from an anisotropic random 
magnetic field 

For a uniform cosmic ray distibution, we will now derive 
the degree of polarization from the anisotropic random 
magnetic field created by the combined effects of com- 
pression and velocity shear in the shock. 

An anisotropic random field in the plane of the sky, 
bt = (b r /,b z r), produces a degree of polarization (Sokoloff 
et al. I 



P = Po 



(&i 3 ) 





Cx 




c 2 



(A. 
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where 

Cx = \{e 2 p + K 2 -l) 

x [(1 — cos i) 2 cos 40—2 cos 20 sin 2 i + (1 + cos i) 2 ] 

+ K 2 cos 20 sin 2 i— \K sin 40 (1 - cos i) 2 , 
C 2 = 2 + {( 2 + K 2 - 1)(1 + sin 2 i cos 2 0) 

— K(K + sin 20) sin 2 i , 

and po ~ 0.7 is the maximum degree of polarization. 

The observed degree of polarization depends on the 
position in the galaxy because of projection effects. As 
expected, the observed degree of anisotropy of the mag- 
netic field components in the sky plane behind the shock 
is independent of azimuth (together with p) for face-on 
view i = 0. The observed anisotropy is maximum when 
the shock occurs at the minor axis = 90° and zero at 
the major axis = when the disc is viewed edge-on, 
i = 90°. 

Appendix B: Parameters of the fitted regular 
magnetic fields 

In Tables IB.ll and IB. 21 we give the parameters of the fit- 
ted regular magnetic field models discussed in Section 0] 
Although a component of the regular field perpendicular 
to the disc plane (B z in Eq. l|B.2|l ) is allowed in the model 
and we searched for fits using this component, a vertical 
field was not required to obtain a good fit in any of the 
rings in either galaxy. Polarization angles at A3. 5 cm and 
A6.2 cm were averaged in sectors with an opening angle 
of 10° in NGC 1097 and 20° in NGC 1365. The greater of 
the standard deviation and the noise within a sector was 
taken as the error in polarization angle. 

The regular magnetic field is modelled as 

B r = R sinpo + Ri sinp! cos(0 — 0i) 

+ R 2 sinp 2 cos(20 - /3 2 ), 
B$ = i?o cospo + Ri cospi cos(0 — (3\) (B-1) 

+ R2 COSP2 COS(20 — fh), 

B z = R z0 + R z \ cos(0 - (3gi) + i? z2 cos(20 - /3 z2 ), 

where Ri is the amplitude of the ?'th mode in units of 
radm -2 , pi is its pitch angle and f3i determines the az- 
imuth where the corresponding non-axisymmetric mode is 
maximum. The magnetic field in each non-axisymmetric 
mode of this model is approximated by a logarithmic spi- 
ral, pi = const, within a given ring. However, the super- 
position of such modes is not restricted to represent a 
logarithmic spiral. For further details of the method, see 
Berkhuijsen et al. I|1997[> and Fletcher et al. I|2L)U4[1 . 

The foreground Faraday rotation due to the magnetic 
field of the Milky Way was neglected in all cases, since 
both galaxies lie far from the Galactic plane; the expected 
(but unknown) contribution of the Milky Way to RM in 
the direction of the two galaxies is about lOradm -2 and 
will not significantly affect the fits. S is the residual of the 
fit and the appropriate \ 2 value, at the 95% confidence 



level, is shown for the number of fit parameters and data 
points. A fit is statistically acceptable if S < \ 2 - The % 2 
values vary from ring to ring (even when the same number 
of fit parameters are used) as some sectors are excluded 
from the model, either because the average measured sig- 
nal is too weak or because the sector represents a clear 
outlier from the global pattern. 
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Table B.l. Parameters of fitted model for NGC 1097, with notation as in Eg. IB. II 



S 



Ring 1.25 < r < 2.5kpc 1.25 < r < 2.5 kpc 2.5 < r < 3.75 kpc 2.5 < r < 3.75 kpc 3.75<r<5kpc 

Azimuth 20° <4>< 190° 200° < cf> < 370° 20° < <j> < 190° 200° < <j> < 370° 0° < <j> < 350° 
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rad m 2 
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deg 
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deg 
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38 
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Table B.2. Parameters of fitted model for NGC 1365, with notation as in Eg. IB. II and ring width 1.75 kpc. 



Ring centre (kpc) 
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Fig. 1. Top row: Total intensity contours and observed E-vectors (£ + 90°) of NGC 1097 at A3. 5 cm (left) and A6.2 cm 
(right) at 6" resolution. The contour intervals are 1, 2, 3, 4, 6, 8, 12, 16, 32, 64, 128 times 40 /LtJy/beam area. The vector 
length is proportional to polarized intensity, 3" length corresponds to 20 ^iJy/beam area. The vector orientations 
are not corrected for Faraday rotation. The background optical image was kindly provided by Halton Arp. Bottom 
row: Polarized intensity contours and observed B-vectors at A3. 5 cm (left) and A6.2 cm (right) at 6" resolution. 
The vector length is proportional to fractional polarization, 3" length corresponds to 20%. The contour intervals are 
1, 2, 3, 4, 6, 8, 12, 16 times 20 /xJy/beam area. The vector orientations are not corrected for Faraday rotation. Here and 
in other maps, the beam size is shown in the bottom right corner of each panel. 
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Fig. 2. Top row: Total intensity contours and observed B-vectors (E + 90°) of the central and southern parts of 
NGC 1097 at A3. 5 cm (left) and A6.2 cm (right) at 10" resolution. The contour intervals are as in Fig. □ with the 
lowest contour level at 80 ^Jy/beam area. The vector length is proportional to polarized intensity, 3" length corresponds 
to 30 /zJy/beam area. The vector orientations are not corrected for Faraday rotation. Bottom row: Polarized intensity 
contours and observed £?-vectors at A3. 5 cm (left) and A6.2 cm (right) at 10" resolution. The contour intervals are as in 
Fig. 2] with the lowest contour level at 20 /iJy/beam area. The vector length is proportional to fractional polarization, 
3" length corresponds to 20%. The vector orientations are not corrected for Faraday rotation. 
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Fig. 3. Left: Total intensity contours and observed 5-vectors (£ + 90°) of NGC 1097 at A6.2 cm at 15" resolution. The 
contour intervals are as in Fig.Q] with the lowest contour level at 150 /iJy/beam area. The vector length is proportional 
to polarized intensity, 3" length corresponds to 30 /iJy/beam area. The vector orientations are not corrected for Faraday 
rotation. The companion galaxy NGC 1097A is visible in the top right corner. Right: Polarized intensity contours and 
observed B-vectors at 15" resolution. The vector length is proportional to fractional polarization, 3" length corresponds 
to 10%. The contour intervals are as in Fig.^ with the basic contour level at 50 /iJy/beam area. The vector orientations 
are not corrected for Faraday rotation. 
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Fig. 4. Top row: Total intensity contours (corrected for primary beam attenuation) and observed B-vectors 
(E + 90°) of NGC 1365 at A3.5 cm (left) and A6.2 cm (right) at 15" resolution. The contour intervals are 
1,2,3,4,6,8,12,16,32,64,128,256,512 times 60 /xJy/beam area. The vector length is proportional to polarized in- 
tensity, 3" length corresponds to 30 /uJy/beam area. The vector orientations are not corrected for Faraday rotation. 
The background optical ESO image was kindly provided by Per Olof Lindblad. Bottom row: Polarized intensity 
contours and observed _B-vectors at A3. 5 cm (left) and A6.2 cm (right) at 15" resolution. The contour intervals are 
1, 2, 3, 4, 6, 8, 12 times 25 /xJy/beam area. The vector length is proportional to fractional polarization, 3" length corre- 
sponds to 10%. The vector orientations are not corrected for Faraday rotation. 
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Fig. 5. Left: Total intensity contours (corrected for primary beam attenuation) and observed i?-vectors (E + 90°) 
of NGC 1365 at A6.2 cm at 25" resolution. The contour intervals are as in Fig. 0] with the lowest contour level at 
80 /ijy/beam area. The vector length is proportional to polarized intensity, 3" length corresponds to 30 /iJy/beam 
area. The vector orientations are not corrected for Faraday rotation. Right: Polarized intensity contours and observed 
-B-vectors at 25" resolution. The contour intervals are as in Fig. with the basic contour level at 40 /iJy/beam area. 
The vector length is proportional to fractional polarization, 3" length corresponds to 10%. The vector orientations are 
not corrected for Faraday rotation. 
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NGC1097 Magnetic Field HPBW=10" - tleprojected 




Fig. 10. Intrinsic orientation of the magnetic field Bp (regular and/or anisotropic random), obtained from polarized 
synchrotron emission in NGC 1097 (left) and NGC 1365 (right), corrected for Faraday rotation and deprojected into 
the galaxies' planes, with the major axis oriented north-south. The vector length is proportional to polarized intensity 
at A3. 5 cm. The resolutions are 10" and 25", respectively; the deprojected beam is shown in the lower right corner. 
North is at top left and top right of the NGC 1097 and NGC 1365 images, respectively. Distances given in axis labels 
are in arcmin. 
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Fig. 11. Top row: Faraday rotation measure (left) and Faraday depolarization (right) of NGC 1097 in grey scales 
between A3. 5 cm and A6.2 cm at 10" resolution, overlaid onto contours of total intensity at A3. 5 cm. Bottom row: 
Faraday rotation measure (left) and Faraday depolarization (right) of NGC 1365 in grey scales at 25" resolution, 
overlaid onto contours of total intensity at A3. 5 cm. 
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Fig. 12. The azimuthal variation of RM in NGC 1097, obtained between A3. 5 cm and A6.2 cm at 10" resolution in the 
rings 1.25 < r < 2.5 kpc (15"-30", left) and 2.5 < r < 3.75 kpc (30"-45", right) where r is the radius from the centre 
of NGC 1097. The error bars show the standard deviation of RM in each sector. Solid lines are RM values obtained 
from the fits to polarization angles described in Sect. 01 with the dotted lines showing the range of RM. The azimuthal 
angle is measured in the galaxy's plane (position angle of the major axis of —45° , inclination of 45° ) counterclockwise 
from the northern major axis. The southern ridge line of polarized intensity is located at 130° azimuthal angle (left) 
and at 160° (right). 
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Fig. 13. The global structure of the regular magnetic field in NGC 1097 (left) and NGC 1365 (right). Both plots are 
face-on views of the modelled galaxy with the major axis shown by a dashed line and the bar axis by a solid line. The 
ring boundaries for NGC 1097 are at r = 1.25, 2.5, and 3.75 kpc , covering the inner and mid-bar, where the half-length 
of the bar is ~ 10 kpc. The NGC 1365 rings are centred at r = 3.5, 5.25, 7.0, 8.75, 10.5, 12.25 and 14 kpc, with a 
half-length of the bar of r ~ 11 kpc. In both figures, the vector lengths are proportional to the modelled Faraday 
rotation measure in a sector and are only shown where RM > 20 rad m~ 2 . The vectors in the outer ring of NGC 1097 
are scaled up by a factor 1.5 compared to those of the inner ring, the vectors for NGC 1365 all have the same scaling. 
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Fig. 14. Top row: Total intensity contours and observed £?-vectors (E + 90°) in the circumnuclear ring of NGC 1097 
at A3.5 cm at 2" (left) and 4" (right) resolutions. The contour levels are 1,2,3,4,6,8,12,16,32,64,128 times the 
basic contour level, which is 150 and 600 /zJy/beam area, in the left- and right-hand panels, respectively. The vector 
length is proportional to polarized intensity, 1" length corresponds to 50 ^Jy/beam area. The vector orientations 
are not corrected for Faraday rotation. The background optical V-band HST image was kindly provided by Aaron 
Barth. Bottom row: Polarized intensity contours and observed B-vectors at A3. 5 cm (left) and A6.2 cm (right) at 4" 
resolution. The vector length is proportional to fractional polarization, 1" length corresponds to 17%. The contour 
levels are 1, 2, 3, 4, 6, 8, 12 times 20 /zJy/beam area. The vector orientations are not corrected for Faraday rotation. 
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Fig. 15. Left: Total intensity contours and observed i?-vectors (E + 90°) in the central region of NGC 1365 at A3. 5 cm 
at 9" resolution. The contour intervals are 1, 2, 3, 4, 6, 8, 12, 16, 32, 64, 128, 256, 512 times 40 /iJy/beam area. The vector 
length is proportional to polarized intensity, 1" length corresponds to 20 /iJy/beam area. The vector orientations are 
not corrected for Faraday rotation. Right: Polarized intensity intensity contours and observed £?-vectors (E + 90° ) of 
NGC 1097 at A3. 5 cm at 9" resolution. The contour intervals are 1,2,3,4,6,8 times 30 /iJy/beam area. The vector 
length is proportional to fractional polarization, 1" length corresponds to 1.7%. - The vector orientations are not 
corrected for Faraday rotation. 
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Fig. 21. The changes in total intensity and degree of polarization arising from a sheared and compressed, initially 
isotropic, random magnetic field as a function of the gas density contrast e p in a shock for the geometry of the galaxy 
NGC 1365. (a): The increase in total synchrotron intensity in the case of constant cosmic ray density, Eq. (|A.6|I . (b): 
The increase in total synchrotron intensity in the case of energy density equipartition between the magnetic field and 
cosmic rays, Eq. (|A.7|) . (c): The degree of polarization in the case of constant cosmic ray density, Eq. I|A.8|) . with 
Po = 0.7. The parameter K = AV r l/(c s d) is the ratio of the turbulent eddy turnover timescale l/c s to the shearing 
timescale d/AV r and is the appropriate measure of the strength of the shearing; high means strong shear. The 
model inclination is i = 40° and the position angle of the shock is <j> = 40°, the parameters of NGC 1365. 



